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ABSTRACT 

Context. Though Arp220 is the closest and by far the most studied ULIRG, a discussion is still ongoing on the main 
power source driving its huge infrared luminosity. 

Aims. To study the molecular composition of Arp220 in order to find chemical fingerprints associated with the main 
heating mechanisms within its nuclear region. 

Methods. We present the first aperture synthesis unbiased spectral line survey toward an extragalactic object. The 
survey covered the 40 GHz frequency range between 202 and 242 GHz of the 1.3 mm atmospheric window. 
Results. We find that 80% of the observed band shows molecular emission, with 73 features identified from 15 molecular 
species and 6 isotopologues. The 13 C isotopic substitutions of HC3N and transitions from H^O, 29 SiO, and CH2CO 
are detected for the first time outside the Galaxy. No hydrogen recombination lines have been detected in the 40 GHz 
window covered. The emission feature at the transition frequency of H31q line is identified to be an HC3N molecular 
line, challenging the previous detections reported at this frequency. Within the broad observed band, we estimate that 
28% of the total measured flux is due to the molecular line contribution, with CO only contributing 9% to the overall 
flux. We present maps of the CO emission at a resolution of 2.9" x 1.9" which, though not enough to resolve the two 
nuclei, recover all the single-dish flux. The 40 GHz spectral scan has been modelled assuming LTE conditions and 
abundances are derived for all identified species. 

Conclusions. The chemical composition of Arp 220 shows no clear evidence of an AGN impact on the molecular 
emission but seems indicative of a purely starburst-heated ISM. The overabundance of H2S and the low isotopic ratios 
observed suggest a chemically enriched environment by consecutive bursts of star formation, with an ongoing burst at 
an early evolutionary stage. The large abundance of water (~ 10 -5 ), derived from the isotopologue H2 8 0, as well as 
the vibrationally excited emission from HC3N and CH3CN are claimed to be evidence of massive star forming regions 
within Arp 220. Moreover, the observations put strong constraints on the compactness of the starburst event in Arp 220. 
We estimate that such emission would require ~ 2 — 8 x 10 6 hot cores, similar to those found in the Sgr B2 region in 
the Galactic center, concentrated within the central 700 pc of Arp 220. 

Key words. Surveys - Galaxies: abundances - Galaxies: active - Galaxies: individual: Arp 220 - Galaxies: ISM - Galaxies: 
starburst 



1. Introduction 

At a redshift of z — 0.018, Arp 220 is the closest ultra- 
luminous infrared galaxy (ULIRG). This galaxy is an ad- 
vanced merger system a s evidenced by the large tidal tails 
observed in the optical (Uoseph fc WrightJll985l ; iKim et all 
2002; iKoda fc Subaru Cosmos Team! l2009tl and th e dou- 
ble nuclei ( separation ~ 0. 98") observed in radio (iNorris 
1988 ). mm (IScoville et "al~lll997t jDownes fc Solomonlll99£ 



Sakamoto et al. 19991) . sub-mm (jSakamoto et al.ll200 



near- IR wavelenghts ( Graha m et al.l 1 19901 : IScoville et al.l 
I1998D . These nuclei are surrounded by two counterrotat- 
ing gas d isks as well as a larger outer disk encom pass- 
ing both (|Sakamoto et al.l 119991: iMundell et al"1l200l[ ). The 
main nuclear powering sources in ULIRGs are thought to 
be starbursts events, active galactic nuclei (AGN), or a 
combination of bot h. Mid-IR surveys of ULIRGs with ISO 
(jGenzel et al.lll998[ ) suggested that 70% -80% of the popu- 
lation are powered by star formation while only 20% — 30% 
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are AGN powered. This result is consistent with that from 
near-IR searches for hidden broad-line regions (BLRs) to- 
wards ULIRGs that lack BLR signatures in the optical 
(|Veilleux et al.lll999D . 

The identification of the main power source becomes 
very elusive in extremely obscured nuclei like those of 
Arp 220. The nuclei in Arp 220 are affe cted by a severe ob- 
scuration at 2.2/xm ([Scoville et al.|[l998t ). Even at 1 mm, the 
dust towards the more massive western nucleus is found to 
be si gnificantly optically thick (r ~ 1, iDownes fc Eckartl 
I2007D . Thus, a number of arguments have been proposed 
in favor of both powering scenarios to explain the large ob- 
served IR luminosity of Arp 220, though none of them have 
been conclusive. 

Among the evidence favoring the AGN-powered sce- 
nario are studies showing hard X-ray emission strongly 
concentrated towards the nuclei, a hard spectrum X-ray 
point source close to the position of the Western nucleus, 
and a softer point s ource towards the Eastern nucleus 
(|Clements et al.l 12002). Additionally, the luminosity ratio 
Lx(2-iokcV)/ Lfir is unusuall y low relative to wha t is ob - 
served in starburst galaxies (jlwasawa et al.l [200ll |2005|) . 
and the 1.98 keV equivalent width of the 6.7 keV line of F e 
is too large to be purely star burst-driven (|Teng et alj|2009j) . 
The high column densities an d therefore obscuration could 
hide a Compton thick AGN (|Sakamoto et al.lll999L [20081: 
IDownes fc Eckartj|2007t ). 

However, the starburst-driven scenario also appears to 
be supported by a number of observations. The detection 
of dozens of radio supernovae (RSNe) in both nuclei of 
Arp 220 is a clear indication of the starburst e vents in this 
galaxy ([Lonsdale et al.ll2006l: IFarra et al.ll2007D . The super- 
nova rate of 4 ± 2yr~ 1 agrees with the derived star for- 
mation rate based on its FIR luminosity (~ 300 M Q yr~ 1 , 
iDopita et al.ll2005T ) , which implies that the starburst traced 
by the detected RSNe is able to pr oduce the observed 
FIR luminosity (jLonsdale et al.l [20061) . However this con- 
clusion is based on assumptions about the truncation of 
the Initial Mass Function (IMF) and that all supernova 
events result in RSNe, whic h requires a dense and com- 
pact starburst environment (jSmith et al.l I1998T ). Towards 
the Eastern nucleus, the agreement between the supernovae 
and the dust surface dens ity suggests a starburst heated 
dust ( Sakamo to et al.ll2008l) . A compari son of the OH mega- 
maser emission ([Lonsdale et al.l |1998|) with the sub-mm 
continuum peak suggests that even if the masers might be 
associated to an AGN , it would not be the m ain contributor 
to the dust heating ('Sakamot o et al.l 1200 8*1. Bright water 
vapor megamaser emission at 183 GHz suggests the pres- 
ence of ~ 10 6 Sgr B2-like hot cores within the central kilo- 
parsec of Arp 220 (jCernicharo et al1l2006bl ). Extended faint 
soft X-ray emission is detected beyond the optical bound- 
aries of the galaxy with bright plumes extending llkpc, 
claimed to be the resu lt of the starburst driven superwinds 
(|McDowell et al.ll2003l) . 

The dense gas phase of the interstellar medium (ISM) 
in Arp220 has been targeted by numerous observations 
of molecular tracers emission in the mm and sub-mm 
wavelenghts. These observations provided additional con- 
straints on the nature of the power source from a ISM 
chemical composition point of view. HNC emission has 
been observe d to be overluminous with respect to HCN 
in Arp 220 (iHuettemeister et a l. 1995; ICernicharo et al] 
20062 lAalto et al.ll2007bh . A similar relative enhancement 



of HNC is only observed in Mrk 231, with a dominant AGN, 
and NGC 4418, with a putative buried AGN, which sup- 
ports the high HNC/HCN ratio to be a chemical indicato r 
of an X-ray dissociation region (XDR, lAalto et al.l 120071)1) . 
The two nuclei of Arp 220 are revealed to show different 
physical condit ions as observed from the high resolution 
maps of HNC (|Aalto et al.l I2009T) . Model calculations at- 
tribute th e observed H.3Q+ emissio n to enhanced X-ray ir- 
radiation ([van der Tak et al~ll2008l ). In the FIR, however, 
the lower level absorption lines in the spectrum of Arp 220 
match those observed in the diffuse clouds in the envelope 
of the Sgr B2 mole cular cloud complex hosting a massive 
star forming event ([Gonzalez- Alfonso et aT1l2004l) . 

Although the molecular line observations in Arp 220 
are still limited to the bri ghter species , the compilation 
of molecular transitions bv iGreve et al.l (|2009l) as well as 
the detection of complex organic species ([Salter et al.| [2008) 
show that this galaxy is one of the brightest molecular emit- 
ters outside the Galaxy together with the starbursts galax- 
ies NGC 253 and M82. This fact turns Arp 220 into a well 
suited candidate for molecular line surveys. 

It is clear that no matter which powering source drives 
the large IR luminosity in Arp 220, it will certainly have an 
imprint on the physical properties and chemical composi- 
tion of the ISM in this galaxy. In this paper we present the 
first unbiased chemical study in the 1.3 mm spectral band 
of Arp 220 with the aim of finding additional physicochem- 
ical clues on the nature of its hidden power source. 



2. Observations and calibration 

Observations were carried out with the Submillimeter 
Array (SMA) in Mauna Kea, Hawaii. A total of 40 GHz 
of the 1.3 mm spectrum of Arp 220 were covered between 
the frequencies of 202 and 242 GHz in the rest frame. This 
range covers the ~ 30% of the 1.3 mm atmospheric window 
and the upper 71% of the SMA 230 GHz receivers nominal 
tunning range. 

The correlator configuration provided an instantaneous 
1.968 GHz bandwidth with a 0.8125 MHz resolution in 
each sideband. Using both sidebands, with a separation of 
10 GHz, yielded 3.936 GHz per frequency setup. To opti- 
mize the frequency coverage, adjacent tuning setups were 
spaced by 1.984 GHz so that the gap between the ob- 
served bands was 16 MHz (~ 21 kms -1 ) except for the fre- 
quencies 212.016 and 231.936 GHz where a gap of 96 MHz 
(~ 127 kms -1 ) was not covered. With the aim of minimiz- 
ing the telescope time required to complete the survey we 
spent only half a track per frequency setup so a total of 
~ 8 GHz could be observed per night. In Fig. Q]we show 
a summary of the frequency tuning setups as well as the 
frequency coverage. 

The phase reference center of the observations was 
aj20oo = 15 h 34 m 57?10 and 5 J2 ooo = +23°30'll'/5 
with frequencies redshiftcd to a velocity of Vlsr = 
5450 km s -1 . Molecular emission towards Arp 220 is ex- 
tended over ~ 2" — 4" in size as obser ved in the CO 
J = 1 — (IDownes fc Solomon! Il998f). J = 2-1 
(ISakamoto et al.l Il999t I Downes fc Eckartl l2007t ). and J = 
3 — 2 (ISakamoto et al] 2008). We used either compact or 
subcompact SMA configurations so that the vast majority 
of this emission would lie within the synthesized beam. 

Gain calibration consisted of 4min observing scans of 
the quasars J1613+342 or J1504+104 at a distance of 13.6° 
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Fig. 1. Summary of the frequency coverage of the line survey between 202 GHz and 242 GHz. The central frequency of 
each band is indicated. The pair of lower (LSB) and upper (USB) sidebands in the observations are indicated by shading 
the same grey scale. The gap between observed bands is 16 MHz but for the two frequencies shown in the center of the 
graphic (at 212.016 and 231.936 GHz) where the coverage gap width is 96 MHz. 

Table 1. Observational details of each frequency setup 



Date 


Frequency 




Config. a 


b 

T225 


Flux 


Gain 


Bandpass 




(BAND) 




(# Ant.) 










2007 May 06 


215.032 GHz 


(USB) 


SC (7) 


0.10 


Titan 


J1613+342 


3c273 




217.016 GHz 


(USB) 


SC (6) 


0.09 


Uranus 


J1613+342 


Uranus 


2007 May 07 


219.000 GHz 


(USB) 


SC (5) 


0.05 


Callisto 


J1613+342 


3c273 


2007 May 08 


220.984 GHz 


(USB) 


SC (5) 


0.05 


Callisto 


J1613+342 


3c273 


2008 Feb 26 


213.048 GHz 


(USB) 


C(7) 


0.16 


Titan 


J1504+104 


3c273 




222.968 GHz 


(LSB) 


C(8) 


0.17 


Callisto 


J 1504+ 104 


J1924-292, Callisto 


2008 Mar 03 


224.952 GHz 


(LSB) 


C(8) 


0.19 


Titan 


J1504+104 


3c273 




226.936 GHz 


(LSB) 


C (8) 


0.20 


Titan 


J1504+104 


J1924-292, Callisto 


2008 Mar 04 


222.968 GHz 


(LSB) 


C (8) 


0.21 


Titan 


J1504+104 


3c273, Titan 


2008 Mar 05 


228.920 GHz 


(LSB) 


C(7) 


0.20 


Titan 


J1504+104 


3c273, Titan 




230.904 GHz 


(LSB) 


C (7) 


0.20 


Callisto 


J1504+104 


3c273, Callisto, J1924-292 



Notes. Array configuration (SC: Subcompact ■ 
225 GHz during observations. 



C: Compact) and number of antennae available. ( ' Average zenith opacity at 



and 14.9° from Arp220, respectively, which were sampled 
every ~ 20 min. The details on the observation of each indi- 
vidual frequency setup are summarized in Table [T] in which 
the observing date, SMA configuration, number of anten- 
nae available, weather conditions, and the calibrators used 
are given. Table [2] presents some parameters of the resulting 
clean maps for each observed frequency. We achieved an av- 
erage rms of 22 mJy/beam on a 6.5 MHz channel width (i.e. 
8 averaged original channels) with synthesized beams rang- 
ing from the smallest (2.9" x 2.2") to the largest (10" x 6.9") 
resulting from compact and subcompact configurations, re- 
spectively. 

Data reduction and calibration were performed using 
the MIR-IDL package and imaging was carried out using 
MIRIAD. 

2.1. Flux calibration and bandpass determination 

Absolute flux density calibration was derived individu- 
ally for each track from observations of Titan, Uranus or 
Callisto as indicated in Table Q] The nominal absolute cal- 
ibration accuracy achieved with the SMA is ~ 15%. As ex- 
plained above, observations were carried out so that there 
were no overlaps among adjacent bands, so the relative flux 
density calibration becomes more important. To account for 
variations of the individual flux density calibrations on each 
of the observed frequency bands we used the flux densities 
measured for the gain calibrators and fitted their emission 
dependance with frequency assuming a power law depen- 
dance as S u oc v a . Given that the observations were all 
carried out within a short period of time (see Table [T|) we 



can assume the flux density of the quasars not to vary sig- 
nificantly (^ 10%). We fitted power laws with a = —1.1 
and -1.6 for J1504 + 104 and J1613 + 342, respectively. 
These power laws are similar to those derived from the up- 
per and lower sidebands on individual observations, which 
supports the assumption of little daily variation on their 
flux densities. Corrections based on deviations of the quasar 
flux densities from the fitted power law were applied to the 
Arp 220 data. These corrections were smaller than 13% in 
absolute flux density in all cases. From these corrections we 
estimate an accuracy in the relative flux density determi- 
nation of ~ 6% across the 40 GHz observed band. 

Bandpass calibration was derived from long integrations 
on 3c273, Uranus, or J1924 — 292. However, the integra- 
tions on Callisto and Titan were included in some cases 
(see Table [lj. Titan shows a strong CO line which appears 
at ~ 235 GHz and thus it was not used for bandpass in 
this setup. The accuracy of the bandpass determination was 
checked on the observed quasars, J1613+342 or J1504+104, 
to look for possible bandpass residuals which might mimic 
line emission. No significant bandpass glitches were found 
that might affect the analysis presented in this paper. We 
used this CO line contamination on Titan to check the con- 
sistency of the relative flux density calibration based on 
the assumption of constant flux density of the quasars. We 
wrongly calibrated the data without masking the channels 
where CO contributes to the total flux measured in Titan 
and subsequently applied the correction derived from the 
power law fitted to J1504 + 104 density fluxes. The result- 
ing corrected flux matched within < 5% the calibrated data 
derived from Titan with the masked CO line. 
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Table 2. Resulting image parameters for each observed frequency 



Frequency 
(GHz) ' 


rrnsic a 
(mjy/beam) 


Synt. Beam 
(B^ aj xBi; in @P.A°) 


s v b 

(mJy) 


T m b/S c 
K/Jy 


Tb/S" 
K/Jy 


203.048 


19.7 


3.9 


X 


3.0 C 


8 +82 


152 


2.54 


7.42 


205.032 


23.5 


8.8 


x 


7.3 € 


8 +76 


135 


0.45 


7.27 


207.016 


21.0 


8.8 


X 


6.2 £ 


8 +76 


143 


0.52 


7.13 


209.000 


22.0 


7.6 


X 


5.0 € 


8 +76 


175 


0.73 


7.00 


210.984 


20.4 


10.0 


X 


6.9 C 


8 +47 


177 


0.40 


6.86 


213.048 


21.2 


3.8 


X 


2.8 € 


8 +82 


145 


2.51 


6.74 


215.032 


24.4 


8.5 


X 


7.0 € 


8 +77 


156 


0.44 


6.61 


217.016 


21.9 


8.4 


x 


6.1 € 


8 +77 


178 


0.50 


6.49 


219.000 


23.3 


7.3 


X 


4.8 € 


8 +76 


255 


0.73 


6.37 


220.984 


21.5 


9.5 


X 


6.7 C 


8 +47 


218 


0.40 


6.26 


222.968 


17.6 


2.9 


X 


2.6 € 


8 -49 


145 


3.22 


6.15 


224.952 


23.3 


3.6 


X 


2.9 € 


8 -49 


165 


2.31 


6.04 


226.936 


20.6 


2.9 


X 


2.6 € 


8 -54 


202 


3.07 


5.93 


228.920 


23.8 


3.2 


X 


2.4 € 


8 -62 


194 


3.02 


5.83 


230.904 


23.9 


2.9 


X 


2.3 € 


8 -61 


503 


3.44 


5.73 


232.968 


18.1 


2.9 


X 


2.4 € 


8 -47 


153 


3.17 


5.63 


234.952 


15.7 


3.7 


X 


2.8 € 


8 -41 


177 


2.16 


5.54 


236.936 


22.7 


2.9 


X 


2.4 € 


8 -50 


204 


3.03 


5.44 


238.920 


25.8 


3.1 


X 


2.4 € 


8 -54 


205 


2.91 


5.35 


240.904 


25.9 


2.9 


X 


2.2 € 


8 -55 


216 


3.30 


5.26 



Notes. ^ a > Calculated for 6.5 MHz channels. ' ' Average flux density measured in the whole ~ 2 GHz band. ' c ' Flux density to syn- 
thesized main beam brightness temperature under the point source assumption. ^' Flux density to source brightness temperature 
assuming a 2" Gaussian emitting region. 



Table 3. Average continuum fluxes in selected line- free 
spectral bands 



Frequency 


Flux 


rms 


# Channels a 


Range 


density 






(GHz) 


(mJy) 


(mJy) 




203.585 - 203.739 


130.3 


28 


31 


203.775 - 204.528 


130.6 


35 


148 


206.420 - 206.963 


132.6 


38 


107 


208.501 - 208.731 


134.1 


58 


46 


212.552 - 212.701 


123.3 


33 


30 


213.643 - 213.894 


132.3 


46 


50 


213.936 - 214.151 


142.0 


31 


43 


214.561 - 214.751 


148.4 


32 


38 


222.465 - 223.146 


136.6 


31 


134 


223.582 - 223.751 


129.0 


29 


34 


224.674 - 225.196 


147.9 


27 


103 


231.409 - 231.849 


153.6 


41 


87 


232.128 - 233.470 


149.6 


29 


263 


234.599 - 235.490 


171.9 


36 


175 


240.063 - 240.662 


178.6 


46 


118 



Notes. {a > Number of averaged 5.1MHz channels to estimate 
the continuum flux. 



3. Continuum emission 

Due to broad spectral features and the prolific molecular 
line emission in Arp220, most of the observing frequency 
setups showed little or no line free frequency ranges to be 
used to subtract the continuum from individual setups. 
Thus, continuum emission was not subtracted in the UV 
data. Fig. [2] shows a composite spectrum of the whole cov- 
ered 1.3 mm band. This figure has been composed with 
the individual spectra extracted from the peak position 



of the unresolved emission at aj20oo = 15 h 34 m 57!22 and 
<W> = +23°30'll'/5. 

Table [2] shows the average flux densities measured in 
each individual observed ~ 2 GHz band, where line emis- 
sion was not subtracted. In order to estimate and subtract 
the continuum emission from the spectra we selected a total 
of 15 narrow frequency bands that we considered as free of 
molecular emission. The positions of these bands and the 
measured continuum flux densities are presented in Tabled 
There we also present the measured rms with respect to 
the mean flux density and the number of 5.1MHz chan- 
nels (~ 6.4 — 7.5 kms -1 ) averaged for each band. These 
values are plotted in Fig. [3] where their associated error 
bars are derived from the measured rms around the mean 
flux density value. The total frequency range covered ap- 
pears to be large enough to show the changes in the con- 
tinuum flux as clearly depicted in the Fig. [3] We fitted a 
power law dependance to the measurements as S v oc v a 
with each measurement weighted as (rms/s/ri) where n is 
the number of averaged channels. A power law index of 
a = 1.3 + 0.2 (dotted line in Fig. [3J is derived by fitting 
all the selected bands in Table EH However, we noticed that 
the flux density measurements at the two highest frequency 
continuum bands significantly disagree (> 10%) with the 
fitted frequency dependance. Such a discrepancy might in- 
dicate that the measured flux densities at the two highest 
frequencies may be overestimated due to a significant line 
contamination or a significant error in the relative flux den- 
sity calibration. The solid line in Fig. [3] represents the fit 
to the data excluding the values at the highest frequencies, 
which results in a slope of a = 1.04 ± 0.02. The difference 
in flux density between both fits is < 5% for the frequency 
range covered in our observations, which is within our cal- 
ibration accuracy (Sect. [2"7lj) . The continuum flux densities 
measured around 212.6, 214.6, and 223.6 GHz also show de- 
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180 - 




210 220 230 240 

Frequency (GHz) 

Fig. 3. Power law fit to the continuum measurements in 
molecular emission free bands across the whole 40 GHz sur- 
veyed (see Table [3]). Fit to all measurements (dotted line) 
result in a dependance as S v oc v . The final fit applied to 
the data (solid line), as derived from the measurements rep- 
resented by filled squares, shows a power law dependance 
as S v oc v . Error bars represent the rms with respect to 
the average within each selected frequency range. See text 
in Sect. for details. 



viations from both fits, but the differences from either of the 
fits are still within < 10%. Finally, we adopted the power 
law v as the frequency dependence for the continuum 
flux densities to subtract it from the spectra. 

The average continuum flux in the 1.3 mm band is 
142 mJy, ranging from 129 mjy to 154 mjy at 203 and 
241 GHz, respectively. 

3.1. Line contribution to total emission flux in broad band 
observations 

The data set presented in this work provides an unique op- 
portunity to accurately estimate the contamination of line 
emission to the continuum flux density measured by broad 
band detectors at millimeter wavelengths. The average to- 
tal flux recovered in the 40 GHz spectral band covered by 
the survey is 195 mjy. As seen in Fig.[5]small subpanel, the 
molecular line emission at these frequencies is dominated 
by the J = 2 — 1 transition of carbon monoxide. However, 
the contribution from other molecular species (see Sect. 
is far from being negligible. Assuming the power law spec- 
trum fitted in the previous section as the continuum level, 
we assess the total line flux of both CO and other species 
by integrating the continuum subtracted spectra. We find 
that the contribution of molecular emission to the overall 
40 GHz measured flux constitutes 28%, from which only a 
9% of the total fl ux is attributed to th e CO J = 2 - 1 
line. The work by lAlbrecht et al.l (|2007f ) estimated an av- 
erage contamination by CO to the 1.3 mm flux density of 
21.3% based on 1.3mm bolometer (~50 GHz bandwidth) 
and CO J — 2 — 1 observations in a sample of 99 galax- 
ies. Although their sample does not include ULIRGs, where 
line contribution appears to be very important, the emis- 
sion of molecules other than CO might be also important 
in some of their sources. In particular, the line contamina- 
tion to the continuum flux density might become significant 



in starbu rsts, where an also prolific molecular emission is 
observed ([Martin et alJl2006ft . 

Our observations show that at short mm and sub-mm 
wavelengths, the flux density correction due to molecular 
line emission is more important than the non-therm al con- 
tribu t ion extrapolated fro m radio measurements (jNorris! 
119881: IScoville et aLI 119971) . For Arp220 the non-thermal 
contribution is estimated to be 7.7mJy (203 GHz) and 
6.7 mJy (241 GHz), equivalent to 6% and 3% of the aver- 
aged total flux at those frequencies, respectively, and well 
below the 28% line contribution. 

3.2. Comparison with previous measurements 

IScoville et all (|1997l) and ISakamoto et al.l (|1999l ) reported 
continuum flux densities of 192 ± 20mJy and 208 mjy, re- 
spectively, in a 1 GHz band centered at 229.4 GHz with the 
OVRO interferometer. Within the same frequency range 
we derived an averaged flux density of 185 ± 3mJy, in 
good agreement with their measurement. However, at these 
frequencies, after subtracting the line contribution of ~ 
39mJy (^ 21% of the total flux density), the effective con- 
tinuum flux yields 146 mjy. The continuum flux density 
would be therefore overstimated by ~ 30% due to the con- 
ta mination from l ine em ission. 

IWoodv et alJ (|l989f ) measured a flux density of 140 ± 
20mJy at a sky frequency of 219.5 GHz also with the 
OVRO interferometer. Luckily this 1 GHz band falls into 
the 223.05 — 224.07GHz rest frequency band where our sur- 
vey shows very little line contamination. Our flux density 
in that range is 146 ± 3mJy, with just an estimated 2% 
of line contribution. In this case, the corrected continuum 
flux density we me asure of 142 mjy agrees with that by 
IWoodv et all (|l989l) . 

The agreement bet ween the flux densit y measured in 
this work and that of IScoville et alJ (|1997l) . with a beam 
~ 6.5 times smaller in area, implies that most of the flux 
density is recovered in both the SMA and OVRO interfero- 
metric observa tions. The discrepan cy between the observed 
flux density bv lWoodv et alJ (Il989h and the single dish flux 
density of 340 ± 80mJv bv lThronson et alJ (|1987l ) is likely 
due to strong line contamination in the frequency tuning 
of the latter (not specified in the publication). Similarly, 
our observations can easily explain the average single dish 
flux of 226 z fc 10 mjy measured in a broad window centered 
at 240 GHz |C arico et aT1ll992f) , where molecular line emis- 
sion such as that of CS J — 5 — 4 and even the J = 3 — 2 
transitions of HCO + and HCN, among many others, will 
strongly contribute to the continuum emission. 

4. Molecular emission 

4.1. Carbon Monoxide 

The subpanel in Fig. [5] shows the prominent emission of 
12 CO J = 2 - 1 with a peak flux of 2.3 Jy. The over- 
all linewidth derived from fitting a single Gaussian to the 
line profile i s 419.3 ± l.Qkms" 1 , in agreement with that 
measured bv lGreve et alJ (|2009h with the JCMT telescope. 
Within the synthesized beam we measure an integrated flux 
of ~ 865Jykms -1 . Fortunately, near the CO line there 
are apparently line-free channels within the observed band, 
and therefore, continuum subtraction could be easily done 
in the UV-plane. Fig. 2] shows the continuum subtracted 
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Fig. 2. Composite of the 1.3 mm frequency scan carried out with the SMA between 202 and 242 GHz. Each individual 
spectrum was extracted from the natural weighted datacubes at the peak of emission position aj20oo = 15 h 34 m 57!22 
and <5j2ooo = +23°30'11'.'5. The calculated continuum level is represented by a straight line (see Sect.[3]for details). The 
CO J — 2 — 1 emission line, labeled in the Figure, dominates the spectral range. This line is shown in a separate box 
in the upper right hand size, where the velocity scale in the X-axis is referred to the rest frequency of the transition. 
The double Gaussian profile fitted to the CO feature (see Sect. 14.11 for details) are shown, both the individual velocity 
component (dotted line) and the total fitted profile (solid line) are shown. 



Table 4. CO J = 2 — 1 fitting parameters 
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Fig. 4. Robust weighted images of the continuum emission 
(grey scale) and CO J = 2 — 1 integrated flux (contours) at 
a resolution of 2.9" x 1.9". Continuum emission has been 
estimated from the 230.904 GHz setup in the ranges not af- 
fected by CO or 13 CS. The stars represent the position of 
the 1.3 mm continuum sources from the high resolution im- 
ages of I Sakamoto et aD (|l999f) . Contours are 3er significant 
levels of 14 mJy beam -1 and 49 Jy km s" 1 for the continuum 
and line flux, respectively. 



integrated CO emission overlaid on top of the line-free 
continuum emission. The total recovered CO flux in the 
map is 1650 ± 50 Jykms" 1 , in good a greement with previ- 
ous single-dish JC MT measurements (|Wiedner et al.112003 : 
iGreve et al.ll2009h . Thus we recover all the single-dish flux 
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3140 ± 30 


5208.4 ± 1.2 
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5461.7 ±2.0 


254.6 ± 1.6 


6.7 



Notes. These parameters are used to constrain the LTE analysis 
of all other molecular species. See text in Sect. 1431 for details. 
Double Gaussian fit with equal linewidth. 



in our maps. This integrated flux is 30% larger than that 
recovered by the higher angular r esolution maps of OVRO 
(1.2" x 1.0" IScoville et al.M l997ft and PdBI (0.7" x 0.5", 
iDownes fc Solomonlll998t ). due to the missing short-spacing 
flux in their maps. Fig.[5]shows the 50 kms" 1 channel maps 
where CO emission is detected over a 800 kms" 1 range, be- 
tween 4950 and 5750 kms -1 . 

The spatial resolution of our observations appears not 
to be enough to disentangle the chemistry of the two nu- 
clei. However, the velocity profile provides us with some 
rough information about the origin of the molecular emis- 
sion. Although the broad emission from both nuclei is sig- 
nificantly blended in velocities, as se en in the high reso- 
lution maps of CO 1 - and 2 - 1 (IScoville et al.l fl997l: 
ISakamoto et al.l [1999), the lower velocity component is 
mostly the emission from the southern region and the west- 
ern nucleus, whereas the higher velocity component mostly 
traces the emission from the northern region and the east- 
ern nucleus. 

We fitted a double Gaussian profile to the spectrum ex- 
tracted from the emission peak position shown in the inset 
in Fig. [5] The fitting assumed both velocity components to 
have equal linewidths. This constraint of double Gaussian 
profile with an equal linewidth had little impact (< 13%) 
in the relative integrated intensities of each component as 
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compared to an unconstrained fitting. The CO parameters 
derived from our double Gaussian fit are shown in Table [4] 
and the fitted profiles are overlaid on top of the CO spec- 
trum in Fig. [2] 



4.2. Line Identification 

We used the rest frequencies provided by the JPL 
dPickett et alJll998h and CDMS (iMtiller et al.ll200lL l2005h 
molecular spectroscopy catalogues to identify each spectral 
feature in the survey. As shown in Fig. [6j a total of 73 in- 
dividual or groups of molecular transitions were identified 
in the 40 GHz frequency band covered. Although an aver- 
age of 1.8 lines/ GHz are detected, given the broad emission 
towards Arp220 of FWHM ~ 330 MHz, as measured from 
CO J = 2 — 1, a significant fraction of the observed 1.3 mm 
band is confusion limited. 

A total of 15 molecular species and 6 isotopical substi- 
tutions were identified and fitted. Among them CH2CO is 
tentatively identified for the first time in the extragalactic 
ISM, as well as the isotopologues 29 SiO, H3, 8 0, and all the 
13 C substitutions of HC3N. Table[S]provides an updated list 
of all the species identified in the extragalactic interstellar 
medium. A total of 46 species and 23 isotopical substitution 
have been identified to date. 



4.3. LTE modelling and fitting of the identified lines 

Similar to what is found in Galactic hot cores, the line con- 
fusion and blending aggravates direct fitting of Gaussian 
profiles to individual spectral features. Thus it becomes 
absolutely necessary to fit synthetic spectra of each identi- 
fied molecular species to the observed spectrum. Synthetic 
spectra of each molecule were calculated using the spec- 
troscopic parameters in the JPL and CDM S catalogues 
(jPickett et alJll998tlMuUer et al.ll200ill2005l) . Assuming lo- 
cal thermodynamic equilibrium (LTE) excitation and opti- 
cally thin emission we can calculate the integrated intensity 
of a given transition for any value of column density (N) 
and excitation tem p eratu re (T ex ) as described in Appendix 
B in iMartm et al.l (|2006l) . Due to both blending and the 
relatively low signal-to-noise ratio for some species we fit- 
ted a double velocity component model with the radial ve- 
locity and linewidth fixed to those derived from 12 CO in 
Sect. 14.11 Therefore, only an automatic exploration in the 
column density and excitation temperature parameters of 
each velocity component was carried out to fit the model 
to the observed spectra. The fitting was performed by min- 
imization of the rms of the residuals after subtracting the 
synthetic spectra from the observations. If the uncertainty 
due to the constraint of double velocity components with 
an equal width, discussed for the 12 CO line fitting, applies 
to other molecular lines, an error in the line integrated area 
of individual components of < 13% may be introduced in 
the synthetic spectral fitting, which is lower than the uncer- 
tainty in the absolute flux density and of the same order as 
the statistical error derived from the Gaussian fit. However, 
even if the integrated intensity of each individual compo- 
nent can be significantly affected, such fitting constraints 
are expected to introduce only a small uncertainty of ~ 1% 
in the overall integrated area of the whole spectral profile 
as estimated from the fit to the 12 CO line. 



As an educated cut-off to the number of lines used to 
generate the synthetic spectra, we used all the listed transi- 
tions with lower level energies E\ < 1000 cm -1 (~ 1440 K) 
in the aforementioned catalogs. Out of the ~ 3000 tran- 
sitions used in the model fitting as shown in Fig. |H1 only 
a limited number contribute significantly to the observed 
spectrum. Table [5] shows the fitted parameters for the sub- 
set of 163 molecular transitions with a peak temperature 
> 25 mK, which corresponds to the average la rms noise 
level at a resolution of ~ 250 kms -1 . These parameters are 
given for both velocity components centered at 5208 and 
5462 kms -1 , which we label as component 1 and 2, respec- 
tively. 

Both Fig.[6]and Table [6] are scaled in br ightness temper- 
ature (7b). As detailed by iDowne 1 lfl989h . the conversion 
factor from flux density (S 1 ) into temperature scale can be 
calculated as 



A 2 
2k' 



Tmb/S — OL ^ s ^ m& 



(1) 



where T m b is the synthesized main beam brightness tem- 
perature and Q s * m b is the solid angle of the convolution 
between the source size and the beam. Usually, for a point 
source we can consider f2 s * m i, ~ &mb- The K/Jy con- 
version factors under this approximation are given in col- 
umn 5 of Table [5] However, in order to account for the 
extended emission of Arp 220 we assumed a source size 
9 S = 2 ", based on the CO J = 2 — 1 maps bv lScoville et al.l 
(fl997h . We can therefore correct for the source beam fill- 
ing factor to calculate the source brightness temperature as 
7b = T m b x fl sirm f>/£l s where £l s is the source solid angle. 
If we use this in Eq. [T] we obtain 



T h /S 



A 2 
2k" 



-fi: 



(2) 



Column 6 in Table [5] presents this conversion factor at each 
of the frequencies covered in this work. We note that this 
conversion factors assume all molecules coexist and there- 
fore no variations in the extent of the emission among 
molecules is considered. 

Table [7] provides a compilation of the source averaged 
physical parameters, namely column densities (N) and ro- 
tational temperatures (T ro t), for each detected species and 
both velocity components. H2 column density has been de- 
rived from the C 18 column density assuming a H2/CO ra- 
tio of 10~ 4 and the isotopical 16 0/ 18 abu ndance ratio of 
150, as derived for the starburst NGC 253 ([Harrison et al.l 
119991) . 

4.3.1. Details on the fitting of individual species 

Here we describe the most relevant details on the fitting of 
particular species. 

CH3OH— The main methanol feature in the observed 
band is the 5 — 4 group at ~ 241.8 GHz. This group of tran- 
sitions together with a number of transitions around 205.7 
and 213.4 GHz allow us to derive rotational temperatures 
of T rot = 15 — 20 K. With this temperature we estimate the 
contribution of all the transitions present in the band. It 
is however difficult to assess the presence of high temper- 
ature methanol due to the blending of most of its higher 
transitions. 
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Fig. 5. CO J = 2 — 1 channel maps in steps of 50kms _1 . First 3 contours are 3<r significant and then 6a up to lOOcr, 
with a — 19mJy. The number at top-left in each panel denotes the LSR velocity in kins -1 . 



CH 3 CCH- With the J = 12 - 11 and J = 13 - 12 
groups of transitions of methyl acetylene we derive rota- 
tional temperatures of T ro t = 16 — 21 K, similar to those 
from methanol. With these parameters we have estimated 
the contribution of the J = 14 — 13 group to the emission 
of the CH 3 CN spectral feature at 239 GHz. 

5*0— Only the 54 — 4 3 and 5s — 44 transitions at 206.1 
and 215.2 GHz are detected unambiguously and not affected 
by blending to other lines. Both have similar upper energy 
levels and therefore result in an unreliable measurement of 
the excitation temperature. The parameters derived from 
these transitions were applied to estimate the 5q — 4s line 
emission at 219.9GHz, observed between the 13 CO and 
C 18 lines. The estimated contribution of this transition 
is not very sensitive to the accuracy of the determined ro- 
tational temperature. 

CH 3 CN- together with HC 3 N, CH 3 CN is the only 
species in the survey showing clear evidence of highly ex- 
cited states. We have detected the emission from CH 3 CN in 
both the ground and the ug = 1 vibrationally excited states. 
We fitted simultaneously the J = 11 — 10 and J = 13 — 12 
i>8 = 0, 1 groups of transitions, as the J = 12 — 11 is par- 
tially blended with 13 CO. The derived excitation tempera- 
tures are ~ 400 — 450 K. Although the uncertainty in the 
excitation temperature is large (20 — 30%), the high energy 
vibrational transitions are confirmed in the three observed 
groups of CH 3 CN transitions. The estimated intensities of 



the J = 12 — 11 group agree well with the observed line pro- 
file intensities and also show clear contribution of the high 
excitation transitions. The groups of transitions at 202.8, 
221.3, and 239.7GHz, identified in Fig. [6j are not shown 
in Table [51 given that these features are due to the contri- 
bution of a few dozens of transitions with expected peak 
intensities at a 0.1 a level. 

H2CO— We find only one line of H2CO unblended, 
namely the 31,2— transition at 225.7 GHz. The 3i 3— 2x,a 
transition at 211.2 GHz appears to be only partially blended 
with the faint 34 SO and the vibrational HC 3 N emission. 
However, these transitions do not have enough dynamic 
range in their upper energy level to get an accurate de- 
termination of the excitation temperature. From the fit to 
these transitions, the contribution of the H2CO transition 
at 218 GHz to the HC 3 N J = 24 - 23 observed line profiles 
was estimated. 

C-C3H2— Although no accurate temperature can be de- 
rived from the multiple transitions of c-C 3 H2, a better fit to 
the observed spectrum is obtained with a high excitation 
temperature of ~ 100 K. The uncertainty in the column 
density derived is less than a factor of two compared to 
that from the fit obtained for a T rot — 20 K. 

HC3N— Vibrationally excited emission of HC 3 N to- 
wards Arp220 has been reported by (Martin-Pintado, in 
prep). Their detection of HC 3 N in the v-j = 1 and, tenta- 
tively, vq = 1 vibrational states are confirmed in all the 
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Table 5. Census of extragalactic molecular species and isotopologues detected 
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Notes. Updated Table from Martin (2009). For references on the first detection of each species see (Martm ct al. 2006) and the 
reported detections of: HD and LiH in absorption at high redshift ijCombes fc Wiklindlll99§ : IVarshalovich et aD l200ll ): H^~ in 
absorption towa rds the ULIRG IRAS 08 572+3915 ijGeballe et alJ 12006): H 3 + in e mission towards the starburst M82 and the 
ULIRG Arp220 (|van der Tak et aLlEfflOl ); CH and NH in absorption towards Arp220 (iGonzalez- Alfonso et all | 2004ri; the tentative 
PNC, N 2 D+ C 33 S, and 13 CH 3 OH towards NGC253 and the Large Magellanic Cloud (iMartm et all l20od: IWane et~aD 120091: 
Martm et ai1l2009bP i; H 2 0+ in emission towards Mrk231 and in absorption towards M82 (|van der Werf et a l. 2010; Wcifi et al.l 
2010] ): OH+ and HF in emission towards Mrk231 ijvan der Werf et al.ll2010l ): and the 13 C isotopologues of HC 3 N, Hf O, 29 SiO, 
and CH2CO being reported in this work. 



four groups of HC3N transitions covered in this survey. 
Moreover, we clearly detect the contribution of the emis- 
sion of HC3N in the vj = 2 state. In order to avoid an 
excessive bias towards very high temperatures in the fit- 
ting, we used only the pure rotational transition and the 
vj = 1, vq = 1, and vj — 2 vibrational states. Subsequently, 
given the high temperatures resulting from the fit, the sig- 
nificant contribution of U4 = 1 and 175 = l/tv = 1 lines 
was calculated from the derived parameters. The fitting to 
this molecule was performed individually for each group of 
ro- vibrational transitions, as the model of all the spectral 
features of HC3N in the survey did not accurately fit the 
observations. Thus individual measurements of the excita- 
tion temperature have been derived for each group. The de- 
rived temperatures averaged for both velocity components 
are ~ 320, 380, and 440 K for the J upper = 23, 25, and 26 
ro-vibrational groups, respectively. The J = 24 — 23 group 
of ro-vibrational transitions are significantly blended with 
the C ls O and H2CO transitions. The excitation tempera- 
ture was fixed to 355 K as estimated for the J up per = 24 
transition from the temperatures derived from the other 
groups. 

HNCO- No unblended transition of HNCO is detected 
so no additional information on this species can be derived 
from this survey. The two brightest HNCO features are 
blended with C 18 and CH3OH. To account for the con- 
tribution of HNCO to the model spectra we considered its 



column density to be a factor of 200 with respect to C 18 
(jMartm et al.lt2009afl . 

C ls O- The J = 2 - 1 transition of C 18 appears in 
one of the most crowded regions of the whole spectrum. Its 
emission is blended with transitions of SO, HNCO, CH3OH, 
and vibrationally excited states of HC3N. We fitted both 
HNCO and C 18 si multaneously, with th eir column den- 
sity ratio fixed f200:l. lMartm et al.ll2""009al) . to the resulting 
spectrum after having subtracted the estimated contribu- 
tion of SO and HC 3 N. 

13 CO- The J = 2 - 1 transition of 13 CO is partially 
blended with the CH3CN J — 12 — 11 group of transitions. 
The transition of 13 CO was fitted after subtracting the es- 
timated CH3CN contribution. 

i?2 8 0— We tentatively identify the spectral feature at 
203.4 GHz as the 3^3 — 22,o transition of H^O. Details on 
the certainty of the identification are given in Sect. 15.2.11 
This is the only species where a single Gaussian fit was 
performed in order to compare with the published result s 
from the main isotopologue, H 2 (jCernicharo et a l. 2006b). 
Fit results are shown in Table [9l 

H 13 CC 2 N- We detected the emission of the three 13 C 
isotopologues of HC 3 N. Both HC 13 CCN and HCC 13 CN 
emit at the same frequencies, so we assumed both species 
to have the same abundance. At least one feature of these 
two isotopologues appears not to be blended to any other 
identified species. H 13 CCCN could not be independently 
fitted due to the blending of the transitions so a similar 
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Fig. 6. Detailed view of the spectral line survey. Spectral resolution of the observed data is smoothed to 20.5 MHz 
(25 — 30 kms -1 across the covered frequency range). The LTE model of the identified molecular species is represented as 
a thick solid curve. The identified molecular features are labeled. 



abundance to the other isotopologues was assumed. For 
the three isotopologues we used a rotational temperature 
of 35 K similar to the one derived from the pure rotational 
transitions of HC3N (see Sect. I5.3|) . 

29 SiO— The spectral feature at 214.3 GHz is tentatively 
identified as the 29 Si isotopical substitution of SiO. This fea- 



ture could have a significant contribution from 13 CH3CN. 
The 30 SiO transition at 211.8 GHz lies in between two ob- 
served frequency setups and therefore we could not iden- 
tify this transition which w as expected to be a factor of 1.5 
fainter than that of 29 SiO (jPenziaslll98l"t ). 
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^CH^CN— The detection of this isotopolgue is unclear 
and therefore is not indicated in Fig. [6] or in the census pre- 
sented in Table [SJ However, from the two groups of tran- 
sitions at ~ 214 and 230 GHz, the latter in the edge of a 
band, we can estimate an upper limit to its abundance. 
We used the same rotational temperatures derived from 
CH3CN for this estimate. The feature at 214.3 GHz can- 
not be accounted for by the emission of this isotopologue. 
Moreover, a lower temperature would change the absolute 
column density derived but does not result in a different 
contribution to the 214.3 GHz line. This limit supports the 
identification of 29 SiO at this frequency. 

CH2CO— The emission of CH2CO could explain part 
of the spectral feature centered at 224.3 GHz and the resid- 
ual emission observed not to be fully fitted by CH3CCH at 
~ 222 GHz. We fixed the excitation temperature to 20 K 
and fitted the rest of lines to get an estimate of the column 
density. Then the contribution to the 222 and 224.3 GHz 
features was calculated. The contribution at 222 GHz per- 
fectly reproduces the observed feature, which supports its 
identification. 



4.3.2. Recombination lines 

Two hydrogen recombination lines are covered by our data 
at rest frequencies of 210.5 GHz (H31a) and 231.9 GHz 
(H30a). 

A detection of H31a towa rds Arp 220 was reported by 
lAnantharamaiah et al.l (|2000h with a peak flux density of 
~ 60 mJy. The width of their reported transition is much 
wider than the transitions they observed at radio frequen- 
cies. Their single dish line profile agrees well with the line 
strength and shape of the emission we observe at 210.5 GHz. 
However, we identify this emission as vibrational transitions 
of HC3N, which we confirm at other frequencies. Indeed, 
lAnantharamaiah et al.l (|2000t) found their observed H31a 
line to be an order of magnitude brighter than predicted 
by their models. Both the line profile and the confirmed 
emission of vibrationally excited HC3N emission makes it 
seem that it was not H31a that they detected. Line contam- 
ination or misidentifications might also be considered for 
the other mm transitions they report at 3mm (H40a and 
H42a), where they also observed the emission to be much 
brighter and broader than expected. This is not surprising 
given that the abundance of molecular transitions reported 
in this work was not known at the time. Our results yield 
that the mm recombination lines towards Arp 220 are much 
fainter that reported, which is in be tter agreement with the 
single ionized component model by lAnantharamaiah et al.l 

mm 

H30a lies right in the middle point between two observ- 
ing setups, in one of the 127 km s -1 gaps in the survey. 
Should this line be bright enough, given the broad emis- 
sion of the recombination line s in Arp 220 of ~ 350 km s -1 
(|Anantharamaiah et al.ll2000D . we might expect to detect 
the edges of the line. However, no signs of emission are ob- 
served in any of the neighboring frequencies which implies 
a non promi nent H30a emission as predi cted by the the 
models from lAnantharamaiah et al. (|2000f) . 



4.3.3. Unidentified features 

The LTE approximation appears to be good enough to fit 
most of the spectral features in the 40 GHz surveyed in 
this study. However there are a number of features which 
are not accurately fitted. There are several significant fea- 
tures which are not yet explained such as the residual emis- 
sion from the CH3CN fit at ~ 221 GHz or that around 
CN at ~ 226 GHz and ~ 227 GHz. Similarly the emission 
at ~ 229.5 GHz is not fitted by the vibrational states of 
HC3N. Additionally, we observe significant residual emis- 
sion between C 34 S and methanol at ~ 241.5 GHz. All these 
residual emission features could be due either to contribu- 
tion from unidentified molecules, or non-LTE emission from 
the identified species, or a departure from the kinematical 
constraints on the two Gaussian components applied in the 
modeling. 

On the other hand, there are other unblended features 
that must be due to unidentified molecular species or transi- 
tions under non-LTE conditions. Such is the case of the fea- 
ture at ~ 212.1GHz in the edge of the 213 GHz frequency 
setup. With a peak intensity of ~ 150 mK, at a ~ 7a level 
above the rms noise, this is likely an unidentified molecular 
species. Even more significant is the double 200 mK peak 
features at 224.0 and 224.4 GHz. The emission from the 
tentatively detected CH 2 CO can partially account for this 
feature but fails to reproduce the broad double-peak emis- 
sion profile. Finally the broad feature at 234.3 GHz, though 
tentative, might be from an unidentified line transition al- 
though we can not rule out the possibility of a residual in 
the continuum subtraction. 

Exploration of a broader band covering more transitions 
and deeper integration would be necessary to accurately 
identify these features. 

5. Discussion 

5.1. Chemical composition of Arp 220 

The column densities of 18 molecular species and iso- 
topologues other than CO derived from the line survey 
in this work allow us to carry out a detailed study of 
chemical composition of the nuclear ISM in Arp 220. In 
Fig. [7] we compare the derived fractional abundances rel- 
ative to H2 in Arp 220 with those in the prototypical star- 
burst NGC253. Molecular hydrogen column densities, as 
detailed in Sect. 14.31 were derived for both sources from 
C 18 Q assuming a 16 O / 18 O=150 as estimated for NGC253 
([Harrison et al.lll999D . For Arp 220 we derive Nn 2 = 1.5 x 
10 23 cm -2 and 2.5 x 10 23 cm -2 for the velocity components 
1 and 2, respectively. The direct comparison between the 
chemical composition of both galaxies show clear overabun- 
dances of H 2 S, HC 3 N, CH 3 CN, and CO+ towards Arp 220. 
Moreover, there is an apparent general overabundance of all 
the detected molecular species in Arp 220 except for H 2 CO. 
By averaging the measured abundances in both galaxies 
for all species but those mentioned with clear differences 
in their relative abundances, we notice that the averaged 
molecular abundances in Arp 220 are a factor of ~ 4.5 and 
~ 3 higher for the 5208 and 5462 kms -1 velocity compo- 
nents with respect the averaged abundance in NGC 253. 
Such a difference might be related to the uncertainty in the 
estimate of the H 2 column density from C 18 in both galax- 
ies. The 16 0/ 18 ratio is subject to a large uncertainty in 
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Table 6. Fitted parameters to the identified transitions Table 6. (Cont.) 



(GHz) 



Molecule 
Transition 



Component 1 
J T B dv T B 
Kkms" 1 mK 



Component 2 
jT B dv T B 



230.538 
207.436 
207.436 
207.438 
207.834 
207.838 
207.838 
241.016 
231.221 
235.380 
235.791 
236.062 
205.791 
213.427 
218.440 
230.027 
239.746 
241.700 
241.767 
241.791 
241.842 
241.879 
241.887 
241.904 
241.905 
205.065 
205.077 
205.081 
222.163 
222.167 
239.248 
239.252 
206.176 
215.221 
219.949 
217.105 
216.710 
211.211 
218.222 
225.698 
213.624 
214.621 
215.974 
216.279 
216.809 
217.532 
217.822 
217.822 
217.940 
218.055 
218.160 
218.449 
218.733 
227.169 
209.230 
209.744 
210.044 
210.527 
210.555 
210.586 
209.573 
209.738 
209.245 
209.362 
210.904 
211.225 
211.226 
211.507 
218.325 
218.861 
219.174 
219.675 
219.707 
219.742 
218.682 
218.854 
218.340 
218.462 
220.070 
220.407 
220.408 
220.700 



CO 2 - 1 

NS 9/2 11/2 - 7/2 9/2 e 
NS 9/2 9/2 -7/2 7/2 c 
NS 9/2 7/2 -7/2 s/2 < 



NS 9/2 11/2 
NS 9/2 9/2 - 



•7/2 9/2 f 
7/2 7/2 f 



NS 9/2 7/2 - 7/2 s/2 f 
C 34 S 5-4 
13 CS 5-4 
CO + 2 3/2 — l 3/2 
CO+ 2 3/2 - l 1/2 
CO+ 2 5/2 - l 3/2 
CH 3 OH li,i - 2 . 2 A+ 
CH3OH li, -Oo.oE 
CH3OH 4 2 2 - 3i 2 E 
CH3OH 3_ 2 2 - 4_i 4 E 
CH3OH 5i, B -4i, 4 A+ 
CH3OH 5 ,5 - 4 0>4 E 
CH3OH 5_i 5 - 4_i 4 E 
CH3OH 5 ,5 - 4 ,4 A+ 
CH3OH 5 2 4 - 4 2 3 A- 
CH3OH 5i, 4 - 4i >3 E 
CH3OH 5 2>3 - 4 2 2 A+ 
CH3OH 5_ 2 ,4 - 4_ 2j3 E 
CH3OH 5 2 , 3 — 4 2 2 E 
CH3CCH 12 2 - 11 2 
CH3CCH 12i - lli 
CH3CCH 12 - llo 
CH3CCH 13i - 12i 
CH3CCH 13 - 12 
CH3CCH 14i - 13i 
CH3CCH 14 - 13 
SO 5 4 - 4 3 
SO 5 5 - 4 4 
SO 5 6 - 4 5 
SiO 5 - 4 
H2S 2 2 — 2i 
H2CO 3i, 3 — 2i 2 
H2CO 3o,3 — 2o,2 
H2CO 3i, 2 — 2i i 
C-C3H2 14 9j6 - 14 8>7 

C-C 3 H2 13 7) 6 — 136,7 

c-C 3 H 2 12 7>6 - 12 6 ; 7 
c-C 3 H2 3 3 ,o — 22,1 

C-C3H2 lis, 6 " H4,7 
C-C3H2 10 5 ,6 - 10 4 , 7 

c-C 3 H 2 61,6 - 5o,s 
c-C 3 H2 60,6 
c-C 3 H 2 5i, 4 
c-C 3 H 2 9 3 ,6 



c-C 3 H 
c-C 3 H2 83,6 

C-C 3 H2 7l,6 
C-C 3 H2 4 3 ,2 

HC 3 N 23 
HC 3 N 23 
HC 3 N 23 
HC 3 N 23 
HC 3 N 23 
HC 3 N 23 
HC 3 N 23 
HC 3 N 23 
HC 3 N 23 
HC 3 N 23 
HC 3 N 23 
HC 3 N 23 
HC 3 N 23 
HC 3 N 23 
HC 3 N 24 
HC3N 24 
HC3N 24 
HC 3 N 24 
HC 3 N 24 
HC 3 N 24 
HC3N 24 
HC3N 24 
HC3N 24 
HC3N 24 
HC 3 N 24 
HC 3 N 24 
HC3N 24 
HC3N 24 



5l,5 
4 2 ,3 
9 2 ,7 
52,4 - 4i, 3 

82,7 

7o,7 
32,1 

22 

22 v 7 = 
22 v 7 = 
22 v 7 = 
22 v 7 = 
22 v 7 = 
22 v 6 = 
22 v 6 = 
22 v 5 = 
22 v 5 = 
22 v s = 
22 v 5 = 
22 v 5 = 

22 v s = 
23 

23 v 7 = 
23 v 7 = 
23 v 7 = 
23 v 7 = 
23 v 7 = 
23 v 6 = 
23 v 6 = 
23 v 5 = 
23 v s = 
23 v 5 = 
23 v 5 = 
23 v 5 = 
23 v 5 = 



ll=lo 

ll=lf 

2i=o 

2l = 2o 
2l = 2f 
ll=lc 

ll=lf 

1/V 7 = 

1/V7 -- 

1/V 7 = 

1/V 7 = 

1/V 7 -- 

1/V 7 = 

ll=lc 
ll=lf 

2i=o 

2l = 2o 
2l=2f 
ll=lo 

ll=lf 

1/V 7 -- 

1/V 7 = 

1/V 7 -- 

1/V 7 -- 

1/V 7 = 

1/V7 ; 



3140 


11600 


17 


.0 


62.8 


13 


.5 


49.7 


10 


.6 


39.2 


17 


.0 


62.6 


13 


,4 


49.6 


10 


,6 


39.1 


100 


.2 


369.5 


31 


.2 


115.1 


14 


.6 


53.8 


26 


.2 


96.8 


25 


.7 


94.8 


15 


_7 


57 9 


14 


.8 


54.7 


14 


4 


53.3 


16 


3 


60 1 


25 


.6 


94.5 


38 


4 


141.6 


9 


.6 


35.5 


7 


.5 


27.6 


17 


2 


63.4 


27 


.2 


100.2 


10 


4 


38.4 


16 


4 


60.7 


9 


.3 


34.4 


31 


.9 


117 6 


30 


x 


111 1 


67 


.0 


247 1 


75 


.9 


280. 1 


115 


.8 


427.1 


61 


.0 


225.2 


51 


.0 


188.4 


59 


5 


219.4 




.9 


25 5 


9 


.0 


33 4 


11 


3 


41.7 


17 


4 


64.1 


13 


4 


49.5 


14 


.9 


55.0 


57 


.0 


210 3 


19 


.0 


70. 1 


37 


.0 


136.3 


15 


2 


56 2 


12 


.3 


45.5 


13 


.6 


50 


g 


.9 


32 8 


22 


_7 


83.8 


231 


.5 


854.2 


89 


.1 


328.5 


89 


1 


328.8 


34 


4 


127 


33 


.8 


124 8 


33 


.8 


124 8 


27 


.5 


101.6 


27 


.6 


101 6 


13 


5 


49.8 


13 


5 


49.8 


13 


5 


50 


13 


x 


48 3 


13 


.1 


48 3 


13 


.6 


50.0 


238 


.7 


880.4 


96 


.2 


354.7 


96 


.3 


355.1 


38 


.9 


143.6 


38 


.3 


141.3 


38 


.3 


141.3 


31 


.5 


116.1 


31 


.5 


116.2 


16 


.0 


58.9 


16 


.0 


58.9 


16 


.0 


59.1 


15 


.5 


57.3 


15 


.5 


57.3 


16 


.0 


59.2 



XI s 


1 


mK 


1820 


c 7 nn 


24 


.7 


Ql 9 


19 


.6 


72.3 


15 


.4 


57.0 


24 


.6 


91.0 


19 


.5 


72.0 


15 


4 


56.9 


107 


.5 


396.5 


72. 


.5 


267.4 


9 


3 


34.4 


46 


,7 


172.2 


84 





309.8 


39 


3 


145.1 


26. 


4 


97.3 


34. 


.3 


126.7 


10 


2 


37.6 


35 


,2 


130.0 


39 




144.2 


55. 




203.2 


76 




280.7 


9 


.5 


35.1 


25. 


.8 


95.1 


9 


.5 


35.1 


17 


2 


63.5 


20. 


.5 


75.5 


6 


.9 


25.4 


19 


.2 


70.8 


27 


.0 


99.7 


13 


.8 


50.8 


19 


4 


71.5 


9. 


4 


34.6 


13 


.2 


48.7 


88 


.1 


325.1 


91 


.6 


337.8 


173. 


.8 


641.0 


128 


3 


473.4 


107 


.6 


396.8 


93 


4 


344.4 


73 


.9 


272.7 


91. 


.6 


337.8 



(GHz) 



Molecule 
Transition 



Component 1 
J T B dv T B 
Kkms -1 mK 



Component 2 

fr B dv t b 

Kkms -1 mK 



15 


9 


58 


.6 


7 


.6 


28 


.1 


9 


4 


34 


.6 


49 


3 


181 


.8 


16 


.4 


60 


.6 


32 


2 


118 


9 


10 


.5 


38 


.8 


10 


.8 


39 


.7 


10 


2 


37 


.6 


7. 


.2 


26 


.6 


20 


.2 


74 


.5 


241 


.6 


891 


.2 


83 


.1 


306 


.5 


83 


2 


306 


.8 


28 


.7 


105 


.9 


28 


.2 


104 





28 


.2 


104 


.0 


22 


4 


82 


.5 


22 


4 


82 


.6 


10 


1 


37 


.3 


10 


.1 


37 


3 


10 


.1 


37 


.4 


9. 


.8 


36 


.0 


9 


8 


36 


.0 


10 


.1 


37 


.4 


141 


.8 


523 


.2 


57. 


.2 


210 


.8 


57 


.2 


211 


.0 


23 


2 


85 


.4 


22 


.8 


84 


.0 


22 


.8 


84 


.0 


18 


.7 


69 


.0 


18 


.7 


69 


.1 


9 


.5 


35 


.0 


9 


.5 


35 


.0 


9 


.5 


35 


.1 


9 


.2 


34 


.0 


9 


.2 


34 


.0 


9 


.5 


35 


.2 



227 .419 
227.977 
228.303 
228.822 
228.859 
228.898 
227.792 
227.971 
227.435 
227.562 
229.236 
229.589 
229.590 
229.891 
236.513 
237.093 
237.432 
237.969 
238.010 
238.054 
236.900 
237.086 
236.184 
236.529 
236.661 
238.401 
238.770 
238.772 
239.082 
219.798 
219.798 
219.798 
241.774 
241.774 
219.560 
220.399 
211.014 
215.839 
202.320 
220.709 
239.096 
208.343 
217.398 
226.454 
235.509 
208.363 
208.363 
208.363 
217.420 
217.420 
217.420 
226.476 
235.532 
202.764 
211.577 
220.390 
229.203 
238.016 
214.385 
202.014 
203.940 
220.178 
222.198 
224.328 
240.187 
226.632 
226.660 
226.664 
226.679 
226.874 
226.875 
226.876 
226.887 
226.892 



HC 3 N 25 - 
HC 3 N 25 - 
HC 3 N 25 - 
HC 3 N 25 - 
HC 3 N 25 - 
HC 3 N 25 - 
HC 3 N 25 - 
HC 3 N 25 - 
HC 3 N 25 - 
HC 3 N 25 - 
HC 3 N 25 - 
HC 3 N 25 - 
HC 3 N 25 - 
HC 3 N 25 - 
HC 3 N 26 
HC 3 N 26 
HC 3 N 26 
HC 3 N 26 
HC 3 N 26 
HC 3 N 26 
HC 3 N 26 
HC 3 N 26 
HC 3 N 26 - 
HC 3 N 26 
HC 3 N 26 
HC 3 N 26 
HC 3 N 26 
HC 3 N 26 
HC 3 N 26 ■ 



24 

24 v 7 : 
24 v 7 : 

24 V 7 : 
24 V 7 : 

24 v 7 : 

24 V 6 : 

24 v 6 : 

24 V 5 : 
24 V 5 : 
24 V 5 : 

24 v 5 : 

24 V 5 : 

24 V 5 : 

25 

25 v 7 : 

25 V 7 : 

25 v 7 : 
25 v 7 : 
25 v 7 : 
25 v 6 ■ 

25 V 6 : 

25 V 4 : 

25 V 5 : 

25 V 5 : 

25 V5 ■ 
25 v 5 - 
25 v 5 : 

25 V 5 : 



HNCO 10io.il 
HNCO lOio.io 
HNCO 10io,9 
HNCO lln ,12 
HNCO llii.n 
C ls O 2-1 
13 CO 2 - 1 
34 SO 5 5 - 4 4 
34 SO 5 6 - 4 5 
CH 3 CN 11 
CH 3 CN 12 3 - 11 3 
CH 3 CN 13 3 - 12 3 



99, 



: ll=le 

■ ll= If 

■ 2i =0 

■ 2i =2o 

' 2l =2 f 

■ ll=le 

■ ll= If 

: l/v 7 -- 
: 1/V 7 -- 
: 1/V 7 -- 
: 1/V 7 - 
: 1/V 7 - 

■ 1/V 7 = 

■ ll=le 

■ ll= If 

■ 2i =0 

■ 2i =2o 
; 2l =2 f 
■■ ll=le 

■ ll= If 
: 1 

- 1/V 7 - 
: 1/V 7 - 
: 1/V 7 - 
1/V 7 -- 
1/V 7 -- 
1/V 7 -- 



10lo.ll 

10io,io 



10 3 



HC 13 CCN 


23 - 


22 


HC 13 CCN 


24 - 


23 


HC 13 CCN 


25 - 


24 


HC 13 CCN 


26 - 


25 


HCC 13 CN 


23 22 


- 22 2 i 


HCC 13 CN 


2323 


- 22 22 


HCC 13 CN 


2324 


- 22 23 


HCC 13 CN 


24 23 


- 23 22 


HCC 13 CN 


24 24 


- 23 23 


HCC 13 CN 


24 25 


— 2324 


HCC 13 CN 


25 - 


24 


HCC 13 CN 


26 - 


25 


H 13 CCCN 


23 - 


22 


H 13 CCCN 


24 - 


23 


H 13 CCCN 


25 - 


24 


H 13 CCCN 


26 - 


25 


H 13 CCCN 


27 - 


26 


29 SiO 5 - 


4 





CH 2 CO 10o,io - 9 ,9 

CH 2 CO 10l, 9 - 9l 8 

CH 2 CO Hi, 11 - 10i,io 

CH 2 CO 11 , 11 - 10o,io 

CH 2 CO Hi, 10 - 10i, 9 

CH 2 CO 12i,i2 - Hi, 11 
CN 2, 



CN 
CN 
CN 
CN 
CN 
CN 



CN 2 



: 3/2- 


-3/2 


~ ll/2,3/2 


; 3/2- 


-5/2 


— ll/2- 


-3/2 


'3/2- 


-1/2 


- ll/2- 


-1/2 


; 3/2- 


-3/2 


— ll/2- 


-1/2 


''5/2- 


-5/2 


- I3/2- 


-3/2 


'■5/2- 


-7/2 


— I3/2- 


-5/2 


'•5/2- 


-3/2 


- I3/2- 


-1/2 


■5/2- 


-3/2 


- I3/2- 


-3/2 


'•5/2- 


-5/2 


- 1.3/2- 


-5/2 



175.2 
71.5 
71.5 
29.3 
28.8 
28.8 
23.8 
23.8 
12.2 
12.2 
12.2 
11.8 
11.8 
12.2 



646.4 
263.7 
263.9 
108.0 
106.4 
106.4 
87.6 
87.7 
44.9 
44.9 
45.0 
43.7 
43.7 
45.1 



123.7 456.3 
55.5 204.8 



55.6 
25.0 
24.7 
24.7 
20.8 
20.8 

11.4 
11.4 
11.4 
11.1 
11.1 
11.4 
7.2 



205.0 
92.3 
91.0 
91.0 
76.6 
76.6 

42.1 
42.1 
42.2 
41.1 
41.1 
42.3 
26.5 



185.0 682.5 
240.6 887.7 



13.5 


49 


.9 


12.5 


46 


.1 


21.1 


78 


.1 


24.4 


89 


.9 


40.1 


148 


.0 


32.5 


119 


.7 


25.8 


95 


3 


20.2 


74 


.7 


12.8 


47 


.1 


13.3 


49 


.2 


13.9 


51 


.4 


10.4 


38 


.2 


10.8 


39 


.8 


11.3 


41 


.5 


25.8 


95 


3 


20.2 


74 


.7 


41.6 


153 


.6 


34.0 


125 


.3 


27.2 


100 


.6 


21.5 


079 


.5 


16.8 


061 


.9 


27.1 


100 


.1 


9.2 


34 





14.0 


51 


.7 


10.4 


38 


.2 


10.0 


36 


.7 


7.0 


25 


.7 


13.8 


50 


.7 


13.9 


51 


.4 


22.1 


81 


.4 


8.3 


30 


.6 



129.9 
59.0 
59.0 
26.9 
26.5 
26.5 
22.4 
22.4 
12.4 
12.4 
12.4 
12.1 
12.1 
12.4 

115.8 
59.0 
59.1 
30.2 
29.9 
29.9 
25.9 
25.9 
8.1 
15.6 
15.6 
15.7 
15.3 
15.3 
15.7 
11.8 
10.6 
9.6 
7.9 
7.2 

302.5 

273.3 
32.5 
60.8 
15.1 
25.5 
29.5 
21.2 
17.1 
13.6 
10.7 

7.0 
7.3 



479.1 
217.5 
217.6 
99.1 
97.7 
97.7 
82.5 
82.5 
45.7 
45.7 
45.8 
44.6 
44.6 
45.9 
427.0 
217.8 
218.1 
111.6 
110.2 
110.2 
95.4 
95.5 
29.9 
57.6 
57.6 
57.8 
56.4 
56.4 
57.9 
43.4 
39.2 
35.5 
29.3 
26.7 
1115.8 
1008.3 
120.0 
224.2 
55.7 
94.3 
109.1 
78.0 
63.1 
50.2 
39.4 

26.0 
27.1 



13. 


.6 


50.2 


10 


1 


39.4 


22 


.0 


81.0 


17. 


9 


66.0 


14 


4 


53.0 


11. 


4 


41.9 


8 


.8 


32.6 


39 


4 


145.0 


25. 


.7 


94.9 


39 


.1 


144.5 


29. 





106.8 


18. 


3 


67.7 


27. 


.8 


102.6 


19. 


.5 


71.8 


11. 


.6 


42.7 


38. 


6 


142.4 


11. 


.5 


42.4 


14. 


.3 


52.8 


39 


.1 


144.3 


62. 





228.6 


23. 


3 


85.8 


7 


.4 


27.3 


7 


.4 


27.1 
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Table 7. Derived LTE physical parameters for each veloc- 
ity component 



Molecule 




N 


T rot 




10 14 


cm- 2 ) 


(K) 




Comp.l 


— Comp.2 


Comp.l — Comp.2 


CO 


16565.8 


- 9527.0 




13 CO 


1327.7 


- 1508.2 




c ls o 


1024.7 


- 1675.3 




H 2 S 


283.2 


- 263.1 




CH 3 CCH 


130.2 


- 70.0 


16 - 21 


CH3OH 


50.4 


- 89.7 


15 - 20 


HC3N a 


29.8 


- 25.8 


36 - 35 


HC 3 N* a 


10.1 


- 8.2 


358 - 376 


CH 2 CO 


9.7 


- 27.1 




SO 


10.1 


- 21.4 


16 - 22 


CH3CN 


7.5 


- 12.1 


392 - 469 


c-C 3 H 2 


9.8 


- 6.1 


117 - 87 


HNCO 


5.1 


- 8.3 




tu 


2.7 


— 8.7 




34 SO 


2.2 


- 7.5 


15 - 24 


CN 


2.5 


- 7.2 




NS 


3.4 


- 5.0 




H 13 CC 2 N b 


5.0 


- 2.7 


35 


H 2 CO 


3.0 


- 4.4 


14 - 16 


c 34 s 


3.6 


- 3.9 




13 cs 


1.1 


- 2.7 




SiO 


1.1 


- 1.9 




29 SiO 


0.4 


- 0.6 




13 CH 3 CN 


< 0.5 


- < 0.8 


390 6 - 470 c 



Notes. HC3N parameters shown for the fitting to the 

ground vibrationally state transitions alone (HC3N) and the 
fit to all observed vibrational transitions (HC3N*). See text 
in Sect. I4XT1 (6) Refers to either H 13 CCCN, HC 13 CCN, or 
HCC 13 CN. (c) These temperatures were fixed to match those 
derived in with the main isotopologues. 



NGC 253 (|Martm et al.ll2010t ) and is likely to be different in 
Arp220. The discrepancy might be even larger if the over- 
produc tion of C 18 Q proposed based on CO isotopologues 
ratios (jMatsushita et al.[|2009[ ) is confirmed. However, as 
discussed in Sect. 15.21 this increase would be marked by a 
significant optical depth in C 18 0. 

In order to avoid the uncertainties derived from the es- 
timate of the H2 column density, we referred the relative 
abundances of the detected species to that of CS. The CS 
molecule trace s the densest gas component within Arp 220 
(n ~ 10 6 cm~ 3 lGreve et alJ [2009). which is likely the source 
of most of the observed chemical complexity. The upper 
panel in Fig. [5] shows a comparison of the abundances rela- 
tive to the isotopologue 13 CS in both Arp 220 and NGC 253. 
The lower panel of Fig. [5] shows the ratio between the abun- 
dances measured in Arp 220 and those in NGC 253. In the 
comparison we assume that 13 CS is only moderately opti- 
cally thick (r < 1) in both galaxies, and that both galax- 
ies have a similar 12 C/ 13 C ratio. Additional uncertainties 
might need to be considered in this comparison if the ex- 
tent of the emission or the considered excitation conditions 
vary among the different species. Thus, only abundance dif- 
ferences above a factor of 2 will be considered as relevant 
between the two sources. Such scatter in the chemical com- 
position are found both in prototypical nearby galaxies and 
Galactic sources as described in Sect. 4.4.1 and 4.4.2 in 
iMartm et all (|2006f ). respectively 

As mentioned in Sect. 14.11 we cannot spatially resolve 
the two nuclei. However, we can infer some information 
from the double Gaussian fit to the line profile, given that 
both nuclei are somewhat separated in velocity. We find 
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Fig. 7. Comparison of the fractional abundances relative 
to H2 derived for the two velocity components in Arp 220 
and those measured in the sta rburst galaxy NGC 253. 
Abundances for NGC 2 53 derived |Martm et all (|2006[ ) and 
H 2 derived from C 18 (|Harrison et al.lll999ft . 



that the abundances in both velocity components, relative 
to either H 2 or 13 CS, are the same within a factor of two. 
But we also find an overabundance of a factor of ~ 4 in 
CH 3 CCH and c-C 3 H 2 , and of a factor of - 2.5 in H 2 S, 
and HC3N in the lower velocity component, presumably 
the western nucleus, 

We find that, with few exceptions, nearly all molec- 
ular species show similar abundances, within a factor of 
2, to those in NGC 253. The main differences, shown by 
both nuclei in Arp 220, are the prominent overabundance 
of H 2 S and the underabundance of H 2 CO in Arp 220 as 
compared with NGC 253. To a lesser extent, we find the 
abundances of HC3N (mostly in the velocity component 1) 
and CO + (for both velocity components) to be a factor ~ 3 
larger in Arp 220. We also find a marginal underabundance 
of molecules such as CH3OH and HNCO, being a factor 
- 2 - 3 lower towards Arp 220 than in NGC 253. Overall 
we find that, with a few exceptions, the chemical compo- 
sition of Arp 220 bears a significant resemblance to that of 
NGC 253. 

5.1.1. Sulfur chemistry in Arp 220 

The molecular abundance that differs the most in Arp 220 
from the starburst galaxy NGC 253 is that of H 2 S. We ob- 
served that all other sulfur bearing molecules in the survey, 
namely SO, NS, and the isotopologues C 34 S and 34 SO, ap- 
pear to have relative abundances similar to those in the 
nuclear region of NGC 253. Our observations show that the 
average relative abundance ratio is (H 2 S/ 13 CS)ai- p 220 ~ 
20 x (H 2 S/ 13 CS)ngc253- However this difference strongly 
depends on the excitation temperature of the species as 
their abundances have been derived from 2 mm transitions 
in NGC 253 and 1.3 mm for Arp 220. In order to evaluate 
whether the observed difference is not due to excitation and 
actually traces an abundance difference, we calculated the 
ratios in both galaxies for different excitation temperatures 
of the gas. Even in the extreme case where the H 2 S would 
trace hot gas (T ex > 50 K), the ratio between the relative 
abundances in both galaxies would drop down to a factor 
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Table 8. Isotopic column density ratios 
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Fig. 8. ( Upper panel) Comparison of the fractional abun- 
dances relative to 13 CS derived for the two velocity compo- 
nents in Arp 220 and those measured in the starburst galaxy 
NGC 253. (Lower panel) Ratio of the relative abundances 
of each velocity component in Arp 220 and in NGC 253. 
Horizontal dotted lines are set to indicate a ratio of factor 
of 2 with the aim of guiding the eye in the comparison. 
Differences below a factor of 2 are not considered represen- 
tative. 



of ~ 3 — 5. It shows that regardless of excitation consid- 
erations, the relative H2S overabundance in Arp 220 with 
respect to NGC 253 would still be significant. 

Though absolute fractional abundances of the sulfur 
bearing species strongly depend on the accuracy of the 
H2 column densities, the relative abundances of these 
spe cies can be reproduc ed by the available chemical mod- 
els (jHatchell et al.l[l998tl . From the comparison to models, 
the overabundance of H2S can have two possible explana- 
tions. The first assumes the grain composition in Arp 220 
to be significantly different from that of NGC 253. The 
H2S frozen in grains, and subsequently injected into gas 
phase, should be one order of magnitude more abundant. 
The other possibility is to assume the grain composition 
to be similar in both objects and con sider time dependant 
sulfur bearing molecules abundances (jHatchell et al.| [l998) 
to be tracing the average age of the star forming molec- 
ular clouds. Though it is somewhat speculative to discuss 
time scales without an accurate measurement of the abso- 
lute abundances, we can try to explain the observed dif- 
ferences in the framework of the evolutionary stage of the 
starburst. In this scenario, the overabundance of H2S would 
imply the hot cores in Arp 220 are, on average, in an earlier 
stage of evolution than NGC 253. The comparison of chem- 
ical models wi th observations of a small sample of massive 
dense cores bv lHerpin et al.l (|2009l ) shows how the CS/H 2 S 
abundance ratio decreases by an order of magnitude as a 
function of the stage of evolution of the clouds. This re- 
sult would also point out the H2S difference observed be- 
tween Arp 220 and NGC 253 as a consequence of evolution- 
ary stage of the starburst. From this observational result, 
the cores in Arp 220, on average, would be in a later stage of 
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evolution than NGC 253, contrary to our previous assertion. 
However, in that case we would not expect such high H2S 
abundance as observed in Arp 220. Therefore, the large rel- 
ative abundance of H2S observed towards Arp 220 appear to 
favour the idea of an earlier stage of starburst evolution in 
this galaxy with respect to that in NGC 253. Nevertheless, 
the abundances of other sulfur-bearing species such as OCS, 
H2CS and SO2 would be ne eded for an accurate comparison 
with the chemical models (lHatchell et al.lll998D. s imilar to 
that performed with NGC 253 (jMartm et al.ll2005l) . 

Regarding the possible AGN-driven origin of the en- 
hanced H2S abundance, there is no evidence in the liter- 
ature to attribute the large abundance of H 2 S to a possi- 
bly hidden AGN. Moreover, recent laboratory experiments 
show the rapid decline of H 2S observed in icy ma ntles when 
exposed to ion irradiation ( Garozzo et al.ll2010f) . 



5.2. Isotopic ratios in Arp 220 

In Table [5] we summarize the isotopologue column density 
ratio derived for each of the species where any isotopic sub- 
stitution has been measured. Ratios are given for each of 
the velocity components as well as the average ratio. In the 
following we will only consider the average ratios, which are 
less affected by the fitting constraints applied for modelling 
the observed spectrum. 

The large 12 CO/ 13 CO ratio (#12/13 > 20) observed in 
luminous mergers as compared to normal starbursts has 
lead to the postulation of an intrinsic differentiation be- 
tween the ISM properties in these objects. An abundance 
difference can be understood if the ISM in mergers is fed 
with unprocessed molecular gas drawn from the external 
region into their nuclear region, or the ISM is enriched 
i n 12 C by nucleos ynthesis in newly formed massive stars 
(|Casoh et al.l fl992). Additionally, these high ratios can oth- 
erwise be explained by different gas surface densities re- 
sulting in differences in opaci ties between lumin ous merg- 
ers and starbursts galaxies (|Aalto et al.l |1995[ ). Towards 
Arp 220, ICasoli et al.l (|1992f) measured a -R 2 ^^ ratio of 
18 ± 5 based on the J = 2 — 1 transition, consistent with 
the limit of R\ ^fi^ > 19 a n d > 22 ratios based on the 
J = 1 — lines (|Aalto et al.lll991t ICasoli et al.lll992[ ). Our 



results give a ratio R 



12/13 



9 towards Arp 220. However, 



this is based on the extracted spectra at the peak of emis- 
sion, and an estimate of the extended emission is required 
to compare with the derived single-dish ratio. As detailed 
in Sect. 14.11 the overall extended CO emission is about a 
factor of 2 larger than that measured towards the central 
position. To estimate the extended 13 CO emission, it would 
require an accurate continuum measurement and subtrac- 
tion as well as precise determination of the emission of the 
other molecules contributing to this crowded spectral band. 
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Thus, with the 13 C0 value at the central peak of emission 
and the integrated extended emission of CO we can con- 
strain the ratio as 9 < < 18. These limits assume 

the extreme situations of 13 CO being as extended as CO or 
compact, respectively. Though the upper limit agrees with 
the avera ge value derived fr om single dish in the literature 
of 18 ± 3 (|Greve et alJl2009h . our limits are consistent with 
the range of the 10 ^ ^1 2/13 & 15 ratio found in normal 
starbursts bv lAalto et all (jl995l) . The 13 CO emission from 
single dish observations may have been underestimated by 
both the baseline subtraction uncertainty due to the lim- 
ited spectral bandwidth and the unknown contamination 
by other species, as well as to a relative calibration issues 
affecting the single dish data. 

The H 12 C 3 N/H 13 C 3 N ratio is - 7, even lower than 
that measured with CO. With a critical density of n cr it ~ 
2 — 5 x 10 6 cm~ 3 , the HC3N emission is likely tracing the 
densest molecular material in Arp 220. Assuming both CO 
and HC3N isotopologue ratios are linked to the carbon iso- 
topic ratio, the similar measured ratio would imply equally 
large opacities affecting the main isotopologues of both 
species. Thus, no reliable estimates of the 12 C/ 13 C isotopic 
ratio can be derived from these species. 

If we assume the SO/ 34 SO ratio of ~ 3 is less affected 
by opacity and is representative of the sulfur 32 S/ 34 S, it 
would result in an extremely low sulfur isotopic ratio, sig- 
nificantly lower than the ratios of 8 — 13 mea sured in the 
starbursts galaxies NGC 253 and NGC4945 (|Wang et al.l 
12004 iMartm et aHl2005l) . This low ratio would be consis- 
tent with the scen ario of 34 S overproduction in supernovae 
(|Chin et al.|[l9 96'). Indeed, t he supernova rate esti mated to- 
wards Arp 220 of 4 ± 2 yr^ 1 (|Lonsdale et al.ll2006l ) is one or 
two orders of magnitude hig her than the rates of 0.05 — 
0.3 yr" 1 towards NGC 253 (|Ulvestad fc Antonuccil Il997t 
iMattila fc Meikld200ll ). However, if we assume this 32 S/ 34 S 
ratio and the measured average ratio of C 34 S/ 13 CS ~ 2, 
we estimate a CS/ 13 CS ratio of ~ 6, similar to the val- 
ues derived from both CO and HC3N. Therefore we sus- 
pect C 34 S and SO to be also affected by significant optical 
thickness effects, resulting in a lower limit to the ratio of 
SO/ 34 SO> 3. 

We measure a 13 CO/C 18 ratio of ~ 1, in agree- 
ment with the ratio derived with the J = 1 — and 
J = 2 — 1 transitions (jGreve et al.l 120091 : iMatsushita et al.l 
2009). Even if we assume an optical depth of 13 CO of 
T13 > 1 (jGreve et al.ll2009D and 12 C/ 13 C iso topic ratio of 
~ 40 — 50, as measured in starburst galaxies (Hcnkel "et al.l 
1993: lHenkel &: MauersbeTg er 1993), we infer a limit to the 
ratio of 16 0/ 18 > 80 - 100. This ratio is lower than that 
derived towards the starburst NGC 253 ( 16 0/ 18 > 150, 
lHarrison et al.llT999! ) , which might imply a large opacity af- 
fecting 13 CO an d even a significant opac ity affecting C 18 
as suggested by IMatsushita et al.l (|2009f) . 

Additionally, we have also detected 29 SiO at a very 
low ratio with respect to the main isotopologue. Though 
the detection remains tentative, we can evaluate the feasi- 
bility of this detection as compared w ith the ratios mea- 
sured in Galactic sources. IWotfll (|1980ll found 28 Si/ 29 Si = 
10 ± 3 toward bright SiO emission regions, namely Sgr A 
and the massive star forming regions Sgr B, Orion A, 
and W51. They observed this ratio to be a factor of 2 
below the terrestrial abundance. Other studies find in- 
tegrated intensity ratios of 5 — 8 and 12 towards the 
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Giant Molecular Clouds c ores W3(H 2 Q) and Sgr B2(M ) 
(jHelmich fc van Dishoeckl Il997t iNummelin et al.l I2000Q . 
and 10 — 20 towards the late carbon star IRC +1026 
(jGroesbeck et al.l ll994; Ka waguchi et al.lfl995Tl . Our detec- 
tion results in a ratio 28 Si/ 29 Si ~ 3 which is on ly a factor o f 
2 below the ratio measured towards Orion A (|Wolfi1ll980l) . 
Even though this result has to be confirmed at other fre- 
quencies and the ratio is subject to a significant error due 
to the low signal-to-noise of the lines involved, the mea- 
sured ratio would tentatively point towards a slight 29 Si 
enrichment in Arp 220. Howe ver this result is ten tative. 

As shown for NGC 253 bv lMartm et al.l (|2010t ). an accu- 
rate measurement of isotopic ratios towards Arp 220 will re- 
quire the observation of doubly substituted isotopologues. 
Our results suggest that the abundances of 34 S and 18 O 
isotope could be similar or even higher than those found in 
local starbursts. A sequence of short and intense starbursts 
(jParra et al.l [2007) or previously enriched gas feeding the 
starburst might be the reason for such enriched ISM to- 
wards Arp 220. However the opacity of the main isotopo- 
logues prevents us from further constraining these isotopic 
ratios. 

5.2.1. H^O: water vapor isotopologue emission 

The spectral feature observed at 203.4 GHz has been identi- 
fie d as the 18 Q isotopologue of the water emission reported 
by ICernicharo et alj (|2006br i. Table M shows the parame- 
ters of the Gaussian profiles fitted to both H2O and H^O. 
The emission of para-H^O 3^3 — 22,o has been detected 
in a number of Galactic hot cores where the emission is 
observed to be significantly blended to the brighter SO2 
12o.i2 — 11 1 ,11 emission ~ 16 MHz below (+24kms _1 ) the 
H^O line (j Jacq et alj [19881: iGensheimer et al.lll996ft . We 
do not find any trace of significant SO2 emission in the fre- 
quency range covered, where a number of SO2 transition, 
brighter than that at 203 GHz, should have been detected. 
Moreover, if the observed feature would be mostly domi- 
nated by SO2 it would appear red-shifted with respect to 
the H 2 velocity. However, the rarer isotopologue emission 
is slightly narrower and blue-shifted (~ — 65kms~ 1 ) with 
respect to t he H2C This is consiste nt with the observed 
H2O profile (|Cernicharo et al1l2006bf ). which is brighter at 
the lower velocities. Also the contribution from methyl for- 
mate and dimethyl ethe r, significant in Galactic hot cores 
([Gensheimer et al1ll996l ). is negligible in Arp 220 based on 
the fitting of these species to the whole 40 GHz band. 

Unlike nearby starburst galaxies, where no H3, 8 obser- 
vations have been reported, the detection of H^O emis- 
sion towards Arp 220 is only possible due to the heavily en- 
riched molecular gas in this galaxy, with very low 16 0/ 18 
as suggested in Sect. 15.21 We derive an integrated flux ratio 
of ^O/H^O ~ 6. This ratio is not representative of the 
16 0/ 18 isotopic ratio given the weak maser origin of the 
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H2O line. Moreover, the ratio between both isotopologues 
might imply that most of the observed H2O emission in 
Arp 220 is indeed a weak maser. 

From the H^O integrated intensity of ~ 70.5 K km s -1 
measured in this work we can estimate a total column den- 
sity Af H i8 = 4.6 x 10 16 cm~ 2 assuming optically thin emis- 
sion and an excitation temperature of T cx = 100 — 200 K. 
Given that we use the C 18 emission as a tracer of the to- 
tal H2 column density, we can reliably derive the fractional 
abundance of water of Xn 2 o — 1.7 x 10 -5 , independent 
of the assumption on the 16 0/ 18 ratio. This value agrees 
well with the fractional abundance ob served towards the 
Sgr B2 region ()Cernicharo et al.ll2006a|) and with an aver- 
age of ~ 10~ 5 towards t he 13 hot cores where H ^O was 
detected in the sample of iGensheimer et alJ (l996). 

We observed a H^O luminosity of 3.9 x 10 7 K km -1 pc 2 
in the centr al 2" (~ 700pc) o f Arp 2 20. The observations 
of H^ 8 by IGensheimer et al.l (|l996l) towards the Sgr B2 
molecular cloud complex shows the emission to be con- 
centrated towards Sgr B2(N) hot core, with a significantly 
lower intensity towards the Sgr B2(M) core. They reported 
an integrated intensity of 17.4 K kms -1 towards Sgr B2(N). 
The emission is unresolved at their 12" (~ 0.5 pc) resolu- 
tion, which is si milar to the 0.4—0.5 pc s ize estimated for the 
northern core (|Lis k. Goldsmith! 1 19901 ). Assuming a range 
in the extent of the emission extent of .5 — 1 pc and us- 
ing a si milar argument to that for H2O in lCernicharo et alJ 
(|2006bh . a total of - 2 - 8 x 10 6 Sgr B2(N)-like cores would 
be required and enclosed in the central 700 pc of Arp 220 to 
explain the observed H^O emissio n. If we consider the lu - 
minosity of Sgr B2(N) of - 10 7 L Q (jde Vicente et al.ll2000ft . 
such number of massive star forming cores alone would eas- 
ily account for the bolometric luminosity of Arp 220. 

Such a concentration of hot cores would require an aver- 
age distance of 3.5— 5.6 pc between cores, assuming they are 
distributed isotropically within a spherical region of 700 pc 
diameter. Assuming them to be distributed within a disk 
with a thickness of ~ 200 pc would yield an average distance 
of 2.4 — 3.7 pc. These rough estimates are still larger than 
the projected distance of ~ 1.6 pc between the hot cores in 
the Sgr B2 molecular complex within the Galactic center 
region. Therefore the estimate of the number of hot cores 
in Arp 220 does not put severe constraints on the compact- 
ness of the star forming regions in this galaxy. However, it 
does imply a widespread intense star formation throughout 
its whole central region. 



5.3. High temperature vibrationally excited molecular gas 

Vibrationally excited HC3N emission probing the hot 
molecular component has been detected, so far, towards two 
extragalactic sources, the luminous inf rared galaxy (LIRG) 
NGC4418 and the ULIRG Arp 220 (iCostariiola fc Aaltol 
I2010L Martin-Pintado, in prep). Bo th galaxies show 
prominent silicate absorption features (jRoche et al.l Il986t 
iSmith et al.l I1989D which can be a sig n of high tempera - 
ture dust heated by hot young stars (jSpoon et al.ll2006h . 
Within the Galaxy, vibrationally excited HC3N emission is 



found towa rds massive (Sgr B2 and Orion Ide Vicente et al 



of the HC3N obse rved transitions are n cr j t > 10 6 cm~ 3 
(jWernli et al. 2007]). Thus, these species are tracing the 
densest and hottest molecular gas in Arp 220. 

In Table [7] two sets of physical parameters are derived 
for HC3N. The first one takes only into account the pure 
rotational transitions of HC3N for which we derive a rota- 
tional temperature of T rot ~ 35 K. The second set of pa- 
rameters with a temperature of T rot ~ 350 K are derived 
with the vibrationally excited transitions only. As observed 
towards the Sgr B2 molecular cloud in the Galactic Center 
region, both HC3N and CH3CN rotational emission arise 
from the same region and trace a m olecular ridge containin g 
the massive star forming regions (|de Vicente et al.lll997l) . 
Vibrational emission, on the other hand, mostly arises from 
the hot cores, though it is not restricted to the core extent 
(|de Vicente et al.l l2000D . In order to populate the vibra- 
tionally excited levels of the observed transitions, HC3N 
and CH3CN require IR photons in the range of 15 — 45 ptm 
and 20 — 30 /im, respectively. While the rotational transi- 
tions arise from the war m envelope around the hot cores 
(40 - 80K, 2 x 10 5 cm- 3 : Ide Vicente et aT1ll997D . the hot 
component will be directly tracing the innermost regions 
around the hot cores w ith temperatures well above 200 K 
(|de Vicente et al.ll2000D . Thus, the excitation of both HC 3 N 
and CH3CN in Arp 220 agrees with its emission originating 
in Sgr B2-like regions. 

6. Conclusions: AGN vs SB driven chemistry 

Up to now, the only available extragalactic molecular 
studies available as compari son templates are the unbi- 
ased study toward NGC 253 (jMartm et al.ll2006h and M 82 
(Aladro et al., in prep), and the targ eted observations to- 
wards NGC 4945 (jWang et al.l 120041) . These studies have 
characterized the chemical complexity of starburst galax- 
ies. The molecular abundances in Arp 220, resulting from 
the detected transitions in the 1.3 mm atmospheric win- 
dow, do resemble those found towards the starburst galaxy 
NGC 253. The only species showing an outstanding abun- 
dance towards Arp 220 is H2S. Although the abundance of 
H2S has only been measured towards NGC 253, this over- 
abundance is likely to be the result of grain disruption and 
subsequent H2S i njection into gas pha se in the early stages 
of star formation (jHatchell et al.lll998l ). This sets the burst 
of star formation in Arp 220 in a very early stage of evo- 
lution and thus, on average, younger than the starburst in 
NGC 253. This comparison suggests that the chemistry in 
Arp 220 is mostly driven by bursts of star formation where 
differences in the observed abundances can be accounted 
for by a difference in the state of evolution of the starbursts 
(jMarti'n et al.ll2006D . 

Unfortunately, no complete chemical templates of the 
ISM in the vicinity of AGNs are available and therefore 
the best tracers in highly X-ray irradiated regions are 
bas ed on chemical models. The results from such mod- 
els (iMeiierink fe Spaanl 120051: IMeiierink et alJl2006l 120071: 
lLoenen et al.l I2008D have motivated a number of obser- 
vational studies focused on the species CO, HCN, HNC, 
HCO+. and CN (iPerez-Beaupuits et alJ 120071: IKrips et all 



2000. I2002H and intermediate (HW2. iMartm-Pintado et al 
2005) star forming regions. Towards Arp 220 we find HC3N 
and CH3CN probing hot molecular gas with derived tem- 
peratures of T cx = 300 — 400 K. The critical densities 



2008; iBaan et alj|2008l: IPerez-Beaupuits et alTl2009D . The 
HCN/HCO+ ratio was used bv (|Aalto et al.ll2007al) to eval- 
uate the presence of an AGN on the LIRG NGC 4418 which, 
like Arp 220, has a heavily ob scured nucleus. Th e HCN and 
HCO + comparative study bv lKrips et al.l (|2008h shows how 
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the HCN/HCO+ ratios in Arp 220 are similar to those mea- 
sured in NGC 1068, the prototypical AGN galaxy. However 
the excitation of HCN differs significantly in both galax- 
ies (jKrips et al J 12008). On the other hand, though over- 
luminous HNC emission in Arp 220 could be attributed 
to either pumping thro ugh mid-IR photon s or the pres- 
ence of an XDR region (|Aalto et al.l l2007bh . the observed 
HN C /HCN ratio suggest a PDR rather than a XDR ori- 
gin (|Baan et al.ll2008l l2010D . We estimate a column density 
ratio CN/HCN~ 0.5, where we have use d the single-dish 
HCN observations by iKrips et al.l (|2008l) and our survey 
measurement of the CN emission. To estimate the HCN 
column density we assumed a source size of 2" similar to 
the one used in this paper. This ratio is lower than the 
range CN/HCN- 1 - 4 that (jPerez-Beaupuits et all [2009) 
claims as indicative of XDR regions. Therefore, the CN 
emission does not appear to be enh anced as observed in 
the c ircumnuclea r disk of NGC 1068 (iGarcia-Burillo et all 
12010(1 . Moreover. iGarci'a-Burillo et al.l (|2010tl found an en- 
hancement of SiO in this region, which does not seem to be 
particularly prominent in Arp 220 as compared to NGC 253. 
Therefore, we do not find any molecular emission that could 
be attributed to X-ray driven chemistry in the nuclear re- 
gion of Arp 220. 

Measured isotopologues ratios appear to point to both 
large opacities affecting the main isotopologues of most ob- 
served species and an enriched molecular material in the nu- 
clear region of Arp 220. ISM enrichment towards Arp 220 is 
similar or even enhanced with respect to that found in other 
starbursts, such as NGC 253 and M82, as a consequence of 
a series of consec utive short and intense starburst events 
(jParra et al.ll20 07). Among the isotopologues, the observed 
H^O luminosity can be accounted for by the emission of 
a few 10 6 hot molecular cores associated with the massive 
star formation within the central 700 pc of Arp 220. Such 
a concentration of cores could be responsible for the whole 
bolometric luminosity of this galaxy, rendering unnecessary 
a significant contribution to the luminosity by a deeply em- 
bedded AGN. Far from the H2O atmospheric absorption, 
the optically thin emission of H^O is proven to be one of 
the best tracers of massive star forming hot cores in highly 
obscured nuclei of galaxies. 

The importance of the star formation within the nu- 
clear region of Arp 220 is further supported by the detec- 
tion of vibrationally excited emission of HC3N and CH3CN, 
with vibrational temperatures > 300 K. Such emission is 
also tracing the molecular component associated with hot 
cores as reported by Martin-Pintado (in prep.). However, 
vibrationally excited HC3N emission with a temperature of 
500 K towards NGC 4418 could be understood as a com- 
pact and deep ly embedded AGN, heatin g up the surround- 
ing material (jCostagliola Sz Aaltd 12010). It is remarkable 
that such vibrationally excited emission has never been re- 
ported towards nearby starbursts. The detection towards 
the ULIRG Arp 220 and the LIRG NGC 4418, might be due 
to the significantly larger contribution of hot core emission 
as a consequence of the higher star formation rates in these 
galaxies. Large opacities affect the pure rotational transi- 
tions of HC^ N CSect. 15.21). Due to the large continuum opac- 
ity at 1mm (|Downes fc Eckartll2007(l . the molecular emis- 
sion we observed cannot arise from the vicinity of an AGN . 
Indeed, the simple models presented lSchleicher et al.l ([2010) 
show the heating sphere of influence of a supermassive AGN 
to be limited to the central ~ 100 pc while the starburst 



heating dominates outside this volume. Therefore, the hot 
gas where the vibrationally excited emission of CH3CN and 
HC3N cannot be tracing the regions around the Compton 
thick AGN but must arise from other regions likely un- 
affected by the AGN radiation and purely heated by the 
starbursting events in Arp 220. 

7. Summary 

In this paper we present the study of a 1.3 mm wavelength 
scan of Arp 220 covering the 40 GHz range between the rest 
frequencies of 202 and 242 GHz. 

Molecular line emission is detected over more than 80% 
of the observed band. Continuum emission, previously over- 
stimated in several works due to line contamination, has 
been estimated from the few line free regions of the spec- 
trum. We derive an average 1.3 mm continuum flux of 
142 mJy, ranging from 129 to 154 mJy across the band. The 
total line contribution to the overall flux measured in the 
40 GHz band is estimated to be ~ 28%, while CO emission 
alone accounts for only ~ 9%. Therefore line contamination 
to thermal continuum is shown to be more critical than the 
non-thermal contribution which only counts for a 3 — 6% of 
the total flux density at these frequencies. 

We present maps of the CO emission detected in the 
Vlsr range between 4950 and 5750 kms -1 . This map likely 
recovers all the single-dish detected flux thanks to the short 
baselines in the SMA sub-compact configuration. 

A total of 73 spectral features are detected, implying a 
detection rate of 1.8 lines/ GHz. 15 molecular species and 
6 isotopologues are identified. The 13 C isotopologues of 
HC 3 N, as well as H^O, 29 SiO, and CH 2 CO, are reported 
for the first time in the extragalactic ISM. The lines in the 
entire 40 GHz spectrum have been fitted to the molecular 
emission of the identified species assuming LTE conditions. 
Of the ~ 3000 transitions used to perform the fit, only 163 
are detected above a 1 a level. 

The molecular composition and abundances derived 
from the modelling favor a purely starburst powering source 
in Arp 220. This conclusion is based on the direct compar- 
ison with the molecular abundances derived for the nearby 
starburst galaxy NGC 253, as well as the detection of the 
water isotopologue, H^O, and vibrationally excited emis- 
sion of HC3N and CH3CN, which are likely tracing the 
hot cores associated with the massive star formation sites 
within the nuclear region of Arp 220. The comparison with 
NGC 253 chemistry shows a prominent overabundance of 
H2S which can be plausibly interpreted as Arp 220 being in 
a very early stage of i ts burst of star forma tion, as predicted 
by chemical models (jHatchell et al.lll998l ). 

Though deep sensitive observations are required to ac- 
curately constrain the isotopic ratios in Arp 220, our mea- 
surements suggest the presence of an isotopic enrichment 
similar to or more prominent than that towards nearby 
starbursts. The large opacities of the main isotopologues 
prevent us from getting further constraints from the ob- 
served isotopologue abundance ratios. Indeed an effect of 
the opacity in HC3N likely similar to that in CO is found. 

From our H^O and C 18 measurements we estimate 
a water abundance of Xh 2 o = 1-7 X 10~ 5 without as- 
sumptions on the 16 0/ 18 isotopic ratio. This abun- 
dance agrees with that derived towards the Galactic ho t 
cores (jGensheimer et all 119961 : ICernicharo et al.l l2006al ). 
Assuming all the H^O emission comes from the massive 
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star forming hot cores within the nuclear region of Arp 220, 
the H^O luminosity of 3.9 x 10 7 K km -1 pc 2 within the cen- 
tral 2" would require a~2-8x 10 6 Sgr B2-like hot cores 
enclosed within its inner 700 pc. The detected vibrationally 
excited transitions of HC3N and CH3CN trace a molecu- 
lar component with temperatures between 300 and 400 K. 
Within the Gala xy, such emission is only obse rved from hot 
molecular cores ()de Vicente et al .1 [20001 1 2 2t) . 

We have not found chemical evidence for the presence 
of a buried AGN in the nuclei of Arp 2 20. Though likely 
prese nt as suggested in various studies (jDownes fc Eckartl 
l2007l and references therein) , the suspected black hole has 
no observable effects in the heating of the nuclear ISM 
in this galaxy. Moreover, the starburst contribution alone 
would be able to explain the huge infrared emission ob- 
served in Arp 220. 
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